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This review covers the important advances in recent years in the physics of thin film ferroelectric 
oxides, the strongest emphasis being on those aspects particular to ferroelectrics in thin film 
form. We introduce the current state of development in the application of ferroelectric thin films 
for electronic devices and discuss the physics relevant for the performance and failure of these 
devices. Following this we cover the enormous progress that has been made in the first principles 
computational approach to understanding ferroelectrics. We then discuss in detail the important 
role that strain plays in determining the properties of epitaxial thin ferroelectric films. Finally, 
we look at the emerging possibilities for nanoscale ferroelectrics, with particular emphasis on 
ferroelectrics in non conventional nanoscale geometries. 
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I. INTRODUCTION 

The aim of this review is to provide an account of the 
progress in the understanding of the physics of ferroelec- 
tric thin film oxides, particularly the physics relevant to 
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B. First-principles investigation of ferroelectric thin filmst 22 



Schottkv barrier formation at mctal-fcrroclcctric iunctionffi fesent and future technology that exploits the charac- 
teristic properties of ferroelectrics. An overview of the 
current state of ferroelectric devices is followed by iden- 
tification and discussion of the key physics issues that 
determine device performance. Since technologically rel- 
evant films for ferroelectric memories are typically thicker 
than 120 nm, characterization and analysis of these prop- 
erties can initially be carried out at comparable length 
scales. However, for a deeper understanding, as well as 
for the investigation of the behavior of ultrathin films 
with thickness on the order of lattice constants, it is ap- 
propriate to re-develop the analysis at the level of atomic 
and electronic structure. Thus, the second half of this re- 
view is devoted to a description of the state of the art in 
first principles theoretical investigations of ferroelectric 
oxide thin films, concluding with a discussion of experi- 
ment and theory of nanoscale ferroelectric systems. 

As a starting point for the discussion, it is helpful to 
have a clear definition of ferroelectricity appropriate to 
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thin films and nanoscale systems. Here we consider a 
ferroelectric to be a pyroelectric material with two or 
more stable states of different nonzero polarization (un- 
like electrets, ferroelectrics have polarization states that 
are thermodynamically stable, not metastable). Further- 
more, it must be possible to switch between the two states 
by the application of a sufficiently strong electric field, 
the threshold field being designated the coercive field. 
This field must be less than the breakdown field of the 
material, or the material is merely pyroelectric and not 
ferroelectric. Because of this switchability of the sponta- 
neous polarization, the relationship between, the electric 
displacement D, and the electric field E is hysteretic. 

For thin film ferroelectrics the high fields that must be 
applied to switch the polarisation state can be achieved 
with low voltages, making them suitable for integrated 
electronics applications. The ability to create high 
density arrays of capacitors based on thin ferroelectric 
films has spawned an industry dedicated to the com- 
mercialization of ferroelectric compute r memories. The 
classi c textbooks on ferroele ctricity (|Lines and Glassl 
ll977|) . l)Fatuzzo and Merl Il967|) . though good, are now 
over twenty years old, and pre-date the shift in emphasis 
from bulk ceramics and single crystals towards thin-film 
ferroelectrics. While much of the physics required to un- 
derstand thin-film ferroelectrics can be developed from 
the understanding of bulk ferroelectrics, there is also be- 
havior specific to thin films that cannot be readily un- 
derstood in this way. This is the focus of the present 
review. 

One of the points that will become clear in the course of 
this review is that a ferroelectric thin film cannot be con- 
sidered in isolation, but rather the measured properties 
reflect the entire system of films, interfaces, electrodes 
and substrates. We also look in detail at the effects of 
strain on ferroelectrics. All ferroelectrics are grown on 
substrates which can impose considerable strains, mean- 
ing that properties of ferroelectric thin films can often 
be considerably different from those of their bulk parent 
material. The electronic properties also have a charac- 
teristic behavior in thin-film form. While bulk ferroelec- 
tric materials are traditionally treated as good insulators, 
as films become thinner it becomes more appropriate to 
treat them as semiconductors with a fairly large bandgap. 
These observations are key to understanding the poten- 
tial and the performance of ferroelectric devices, and to 
understanding why they fail when they do. 

In parallel with the technological developments in the 
field, the power of computational electronic structure the- 
ory has increased dramatically, giving us new ways of un- 
derstanding ferroelectricity. Over the last fifteen years, 
more and more complex systems can be simulated with 
more accuracy; and as the length scales of experimental 
systems decrease, there is now an overlap in size between 
the thinnest epitaxial films and the simulated systems. 
It is therefore an appropriate and exciting time to review 
this work, and to make connections between it and the 
problems considered by experimentalists and engineers. 



Finally we look at some issues and ideas in nano-scale 
ferroelectrics, with particular emphasis on new geome- 
tries for ferroelectric materials on the nanoscale such as 
ferroelectric nanotubes and self-patterned arrays of fer- 
roelectric nano-crystals. 

We do not attempt to cover some of the issues which 
are of great importance but instead refer readers to re- 
views by other authors. Some of the more important ap- 
plications for ferroelectrics make use of their piezoelec- 
tric properties, for example in actu ators and micr osen- 
sors; this topic has been reviewed by Muralt (2000). Re- 
laxor ferroelectrics in which ferroelectric ordering occurs 
through the interaction of polar nanodomains induced by 
substitution are also of great interest for a num ber of ap- 
plicat ions and have recently been reviewed by ISamaral 
(200jJ. 



II. FERROELECTRIC ELECTRONIC DEVICES 
A. Ferroelectric Memories 

The idea that electronic information can be stored in 
the electrical polarization state of a ferroelectric material 
is a fairly obvious one; however it's realization is not so 
straightforward. The initial barrier to the development 
of ferroelectric memories was the necessity to make them 
extremely thin films, because the coercive voltage of fer- 
roelectric materials is typically of the order of several 
kV/cm, requiring sub-micron thick films to make devices 
that work on the voltage scale required for computing 
(all Si devices work at < 5V). With today's deposition 
techniques this is no longer a problem, and now high- 
density arrays of non-volatile ferroelectric memories are 
commercially available. However, reliability remains a 
key issue. The lack of good device models means that 
design of ferroelectric memories is expensive and that 
it is difficult to be able to guarantee that a device will 
still operate ten years into the future. Because compet- 
ing non-volatile memory technologies exist, ferroelectric 
memories can succeed only if these issues are resolved. 

A ferroelectric capacitor, while capable of storing in- 
formation, is not sufficient for making a non- volatile com- 
puter memory. A pass-gate transistor is required so that 
a voltage above the coercive voltage is only applied to 
the capacitor when a voltage is applied to both the word 
and bit line (this is how one cell is selected from an array 
of memories). The current (measured through a small 
load resistor in series with the capacitor) is compared to 
that from a reference cell that is poled in a definite di- 
rection. If the capacitor being read is in a different state 
the difference in current will be quite large (the displace- 
ment current associated with switching accounts for the 
difference). If the capacitor does not switch because it 
is already in the reference state, the difference in cur- 
rent between the capacitor being read and the reference 
capacitor is zero. 

Most memories use either a 1T-1C design (1 
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FIG. 1 a)lT-lC memory design. When a voltage is applied 
to both the word and bit line, the memory cell is addressed. 
Shown also is the voltage applied to the capacitor and the cur- 
rent output, depending on whether a one or a zero is stored. 
The current for the zero state is pure leakage current and by 
comparison to a reference capacitor can be removed, (b) A 
2T-2C memory cell in which the reference capacitor is part of 
the memory cell 



Transistor- 1 Capacitor) or a 2T-2C (2 Transistor-2 Ca- 
pacitor) design (Fig. . The important difference is that 
the 1T-1C design uses a single reference cell for the en- 
tire memory for measuring the state of each bit, whereas 
in the 2T-2C there is a reference cell per bit. A 1T-1C 
design is much more space-effective than a 2T-2C design, 
but has some significant problems, most significantly that 
the reference capacitor will fatigue much faster than the 
other capacitors, and so failure of the device occurs more 
quickly. In the 2T-2C design the reference capacitor in 
each cell fatigues at the same rate as its correspond- 
ing storage capacitor, leading to better device life. A 
problem with these designs is that the read operation is 
destructive, so every time a bit is read it needs to be 
written again. A ferroelectric field effect transistor, in 
which a ferroelectric is used in place of the metal gate on 
a field effect transistor, would both decrease the size of 
the memory cell and provide a non-destructive read out; 
however, no commercial product has yet been developed. 
Current efforts seem to run into serious problems with 
data retention. 

An example of a real commercially available memory is 
the Samsung lead zirconate titanate based 4 Mbit 1T-1C 
ferroelectric memory. The SEM cross-section (Fig. 0) 
of the device gives some indication of the complexity of 
design involved in a real ferroelectric memory. 

Lead zirconate titanate (PZT) has long been the 
leading material considered for ferroelectric memories, 
though strontium bismuth tantalate (SBT), a layered 
perovskite, is also a popular choice due to its superior 
fatigue resistance and the fact that it is lead free(Fig. [3Jl. 
However it requires higher temperature processing, which 
creates significant integration problems. Some recent 
progress has been made in optimizing precursors. Un- 
til recently the precursors for Sr, Bi, and Ta/Nb did not 
function optimally in the same temperature range, but 
last year Inorgtech developed Bi(mmp)3 - a 2-methoxy- 
2-propanol propoxide that improves reaction and lowers 
the processing temperature for SBT, its traditional main 
disadvantage compared to PZT. This material also sat- 
urates the bismuth coordination number at 6. Recently 




FIG. 2 Cross-sectional SEM image of the Samsung 4Mbit 1T- 
1C 3 metal FRAM 



several other layered perovskites, for example bismuth 
titanate, have also been considered. 
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FIG. 3 (a)ABOz cubic perovskite structure, (b) Strontium 
Bismuth Tantalate (layered perovskite structure) 

As well as their applications as FRAM's, ferroelectric 
materials have potential use in DRAM's because of their 
high dielectric constant in the vicinity of the ferroelec- 
tric phase transition, a topic which has been reviewed by 
iKingon. Maria and Streifferl l)200fi) . Barium strontium 
titanate (BST) is one of the leading materials in this 
respect since by varying composition a transition tem- 
perature just below room temperature can be achieved, 
leading to a high dielectric constant over the operating 
temperature range. 



B. Future prospects for non-volatile Ferroelectric Memories 

There are two basic kinds of ferroelectric random ac- 
cess memories in production today: (1) The free-standing 
RAMs and (2) fully embedded devices (usually a CPU, 
which may be a CMOS EEPROM (complementary metal 
oxide semiconductor electrically erasable programmable 
read-only memory, the current generation widely used 
non- volatile memory technology), plus an FRAM (fer- 
roelectric random access memory), and an 8- bit micro- 
processor). The former have reached 4 Mbit both at 
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Samsung (using PZT) and Matsushita (using SBT). The 
Samsung device is not yet, as far as the authors know, in 
commercial production for real products, but the NEC 
FRAM is going into full-scale production this year in 
Toyama (near Kanazawa). Fujitsu clearly leads in the 
actual commercial use of its embedded FRAMs. The Fu- 
jitsu embedded FRAM is that used in the SONY Playsta- 
tion 2. It consists of 64 Mbit of EEPROM plus 8 kbit of 
RAM, 128 kbit ROM, and a 32-kbit FRAM plus security 
circuit. The device is manufactured with a 0.5-micron 
CMOS process. The capacitor is 1.6 x 1.9 microns and 
the cell size is either 27.3 square microns for the 2T-2C 
design or 12.5 square microns for the 1T-1C. 

The leading competing technologies in the long term 
for non volatile computer memories are FRAM (ferro- 
electric random access memories) and MRAM (magnetic 
random access memories). These are supposed to replace 
EEPROMs (electrically erasable programable read-only 
memories) and "Flash" memories in devices such as dig- 
ital cameras. Flash, though proving highly commercially 
successful at the moment, is not a long term technology, 
suffering from poor long term endurance and scalability. 
It will be difficult for Flash to operate as the silicon logic 
levels decrease from 5V at present to 3.3V, 1.1V, and 
0.5V in the near future. The main problem for ferro- 
electrics is the destructive read operation, which means 
that each read operation must be accompanied by a write 
operation leading to faster degradation of the device. The 
operation principle of MRAMs is that the tunneling cur- 
rent through a thin layer sandwiched between two fer- 
romagnetic layers is different depending on whether the 
ferromagnetic layers have their magnetization parallel or 
anti-parallel to each other. The information stored in 
MRAMs can thus be read non-destructively, but their 
write operation requires high power which could be ex- 
tremely undesirable in high density applications. We 
present a summary of the current state of development 
in terms of design rule and speed of the two technologies 
in the following table. 



Company 


Design Rule 


Speed 




(Feature Size) 


(Access Time) 


MRAMs 






NEC/Toshiba 


1Mb 




IBM 


16 Mb 




Matsushita 


4Mb 




Sony 


8kb 0.18 microns 




Cypress 


256 kb 


70 ns 


State of the art 


16 Mb 0.09 microns 


25 ns 


FRAMs 






Fujitsu 


32 kb 


100 ns 


Samsung 


32 Mb 0.18 microns 


60 ns 


Matsushita 


4 Mb 


60 ns 


Laboratory 




800 ps 



Some clarification of the numbers in this table is re- 
quired. The size of the Fujitsu FRAM memory may seem 



small but it is for an actual commercial device in large 
scale use (in every Playstation 2), whereas the others 
are figures from internal sampling of unreleased devices 
that have not been commercialized. No MRAMS exist 
in any commercial device, giving FRAMs a substantial 
edge in this regard. The most recent commercial FRAM 
product actually shipped is a large-cell-area six-transistor 
four-capacitor (6T-4C) memory for smart credit cards 
and radio frequency identificatio n tags (RF-ID) and fea- 
tures non-destructive read out l|Masui et all 12003ft . A 
total of 200 million ferroelectric memories of all types 
have been sold industry-wide. The Sony MRAM, though 
small, has sub-micron design rules, meaning that in prin- 
ciple a working device could be scaled up to Mb size. 

Partly in recognition of the fact that are distinct 
advantages for both ferroelectrics and ferromagnets, 
there has been a recent flurry of activity in the field 
of multiferroics, i.e. materials that display both ferro- 
electric and magnetic ordering, the hope being that one 
could develop a material with a strong enough coupling 
between the two kinds of ordering to realize a device 
that can be written electrically and read magnetically. 
In general multi-ferroic materials are somewhat rare, 
and certainly the conventional ferroelectrics like PbTi03 
and BaTiOs will not display any magnetic behaviour 
as the Ti-0 hybridisation required to stabilize the 
ferroelcctricity in these compounds will be inhibited by 
the partially fil led d-orbita ls that would be required 
for mag netism ijHill I |2000|) . However there are other 
mechanisms for ferroelectricity and in materials where 
ferroelecticity and magnetism co-exist there can be 
coupling between the two. For example, in BaMnF4 
the ferroelectricity is actually responsible for changing 
the antiferro magnetic ordering to a weak canted fer- 
romagnetism llFox et all . ^80) In addition, the large 
magnetoelectric coupling in these materials causes 
large dielectric anomalies at the Neel Temp erature 
and a t the in - plane spin ordering temperature l|Scott I 
ll977h . l|ScottL Il979l) . More recent theoretical and 
experimental effort s ha ve fo c used on BiMnO^BiFeO^ 
llSeshadri and HillT (200l[) . llMoreira de Santos et al. , 
2002 1) . tWang et all 1200.1) and YMnO^ l|Van Aken et al . 
2004 . (|Fiebig et alJ . l2002|) . 



C. Ferroelectric FET's 

It has been known for some time that replacing the 
metal gate in a field effect transistor (FET) by a fer- 
roelectric could produce a device with non-destructive 
read-out (NDRO), in which the polarization of the 
gate (+ for "1" and - for "0") could be sensed sim- 
ply by monitoring the source-drain current magnitude. 
Thus such a device requires no reset operation after 
each READ and will experience very little fatigue in 
a normal frequent-read, occasional-write usage. The 
early ferroelectric FE Ts utilized g a tes of lithium nio- 
bate (Rice Univ.) ijRabson et all Il995|) or BaMgF 4 
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(Wes ti nghouse) ijSinharov et all ll99l[) . l)Sinharov et all 
19921. iSinharov et all 11993). An example of a ferroelec- 
tric FET device as fabricated bv lMathews et al\ l|l997h 
is shown in Fig. 0] 
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FIG. 4 \Mathews et alJ ^,997) Schematic diagram of an all- 
perovskite ferroelectric FET and measurement circuit 

The optimum parameters for such a ferroelectric gate 
material are extremely different from those for pass-gate 
switched capacitor arrays; in particular, the latter re- 
quire a remanent polarization ca. 10 /kC/cul 2 , whereas 
the ferroelectric-gated FETs can function well with 50x 
less (0.2 /iC/cm 2 ). However, the switched capacitor array 
(FRAM) is very tolerant of surface traps in the ferroelec- 
tric (which may be ca. 10 20 cm~ 3 in the interface region 
near the electrode), since the ferroelectric makes contact 
only with a metal (or metal-oxide) electrode. By com- 
parison, the ferroelectric gate in an FET contacts the 
Si substrate directly (MOSFET channel - metal-oxide- 
semiconductor field-effect-transistor channel). Thus it 
must be buffered from the Si to prevent charge injection. 
Unfortunately, if a thin buffer layer of a low-dielectric 
material such as SiC>2 is used, most of the applied volt- 
age will drop across the buffer layer and not the ferro- 
electric gate, making it impossible to switch the gate. 
As a result, much of the ferroelectric FET research has 
employed buffer layers with relatively high dielectric con- 
stants, or else rather thick buffer layer s, for exampl e, the 
first BaMnF 4 FET made at Symetrix (jScotd. 11998^ used 
a buffer layer of ca. 40 nm of Si02- Subsequent studies 
often used PZT. ljKalkur et ait 119941 although the large 
remanent polarization in this case (ca. 40 ^C/cm 2 ) is 
actually undesirable for a ferroelectric FET gat e. 

As pointed out bv lYoon and Ishiwaral (|200lt) . the de- 
polarization field in a ferroelectric gate is inevitably gen- 
erated when the gate is grounded, and this makes it 
very difficult to obtain > 10 year data retention in an 
FE-FET. Their solution is to utilize a 1T-2C capaci- 
tor geometry in which this depolarization field is sup- 
pressed by poling the two capacitors in opposite direc- 
tions. With this scheme Ishiwara and his colleagues 
achieved an on/off source-drain current ratio of >1000 
for a 150 nm thick SBT film in a 5 x 50 micron MOS- 
FET channel, with Pt electrodes on the SBT capacitor. 

Note that the direct contact of the ferroelectric onto 
Si produces a semiconductor junction that is quite differ- 
ent from the metal-dielectric interface discussed above. 



The Schot tky barrier heights for th i s case have been cal- 
culated bv lPeacock and Robertson! {2002). The electron 
screening length in the Si will be much greater than in the 
case of metal electrodes; in particular this will increase 
the minimum ferroelectric film thickness required to sta- 
bilize the device against depolarization instabilities. Al- 
though this point was first emphasized by Batra and Sil- 
verman (1973), it has been neglected in the more recent 
context of ferroelectric FETs. In our opinion, this depo- 
larization instability for thin ferroelectric gates on FETs 
is a significant source of the observed retention failure 
in the devices but has not yet been explicitly modeled. 
If we are correct, the retention problem in ferroelectric 
FETs could be minimized by making the ferroelectric 
gates thicker and the Si contacts more conducting (e.g., 
p+ rather than p) . See Scott (Microelectron. Eng. 2005, 
in press) for a full discussion of all-perovskite FETs. 

Table I lists a number of the most promising gate ma- 
terials under recent study, together with the buffer lay- 
ers employed in each case. Studies of the I(V) charac- 
teristics of such ferroel ectric FETs have been given by 
iMacleod and Hoi (l200lt) and a disturb-fre e programming 
scheme described bv lUllman et al\ l)200lt) . 
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Beyond its use in modulating the current in a semicon- 
ductor channel the ferroelectric field effect can also be 
used to modif y the properties of more exot i c corr elated 
oxide systems. ijAhn.Triscone and Mannhardl2003|) 



D. Replacement of gate oxides in DRAMs 

At present there are three basic approaches to solving 
the problem of SiC>2 gate oxide replacement for DRAMs: 
The first is to use a high-dielectric ("high-k") mate- 
rial such as SrTiOs (k = 300 is the dielectric constant; 
e = k — 1 is the permittivity; for k >> 1 the terms 
are nearly interchangeable) deposited by some form of 
epitaxial growth. This is the technique employed at Mo- 
torola, but the view elsewhere is that it is too expen- 
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sive to become industry process-worthy. The second ap- 
proach is to use a material of moderate k (of order 20), 
with Hf02 favored but ZrC>2 also a choice. Hafnium 
oxide is satisfactory in most respects but has the sur- 
prising disadvantage that it often degrades the n-channcl 
mobility catastrophically (by as much as xl0,000). Re- 
cently ST Microelectronics decided to use SrTiC>3 but 
with MOCVD deposition from Aixtron, thus combining 
high dielectric constant and cheaper processing. 

The specific high-k integration problems are four: (1) 
depletion effects in the polysilicon gate; (2) interface 
states; (3) strain effects; and (4) etching difficulties (HTO2 
is hard to wet-etch) . The use of a poly-Si gate instead of 
a metal gate produces grain boundary stress in the poly, 
with resultant poor conductivity. This mobility degrada- 
tion is only partly understood. The general view is that 
a stable amorphous HfC>2 would be a good strain-free so- 
lution. Note that Hf02 normally crystall izes into two or 
three phases, one of which is monoclinic ijMorrison et all 
2003). Hurley in Cork has been experimenting with 
a liquid injection system that resembles Isobe's earlier 
SONY device for deposition of viscous precursors with 
flash evaporation at the target. 



III. FERROELECTRIC THIN FILM DEVICE PHYSICS 

We now turn to some of the physics questions which 
are relevant to ferroelectric thin film capacitors. 



A. Switching 

In the ferroelectric phase ferroelectric materials form 
domains where the polarisation is all aligned in the same 
direction, in an effort to minimise energy. When a field 
is applied the ferroelectric switches by the nucleation of 
domains and the movement of domain walls and not by 
the spontaneous reorientation of all of the polarisation 
in a domain at once. In contrast to ferromagnets where 
switching usually occurs by the sideways movement of 
existing domain walls, ferroelectrics typically switch by 
the generation of many new reverse domains at particular 
nucleation sites, which are not random; i.e., nucleation 
is inhomogeneous. The initial stage is nucleation of op- 
posite domains at the electrode, followed by fast forward 
propagation of domains across the film, and then slower 
widening of the domains(Fig. |SJ. In perovskite oxides 
the final stage of the switching is usually much slower 
than t he other two stages, as first established by iMerd 
In other materials nucleation can be the slowest 
(rate-limiting) step. 



1. Ishibashi-Orihara Model 
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FIG. 5 The three phases of domain reversal, I. Nucleation 
(fast) II. Forward growth (fast) III. Sideways growth (slow) 



model l|Orihara et all Il994h based on Kolomogorov- 
Avrami growth kinetics. In this model one considers 
a nucleus formed at time t' and then a domain prop- 
agating outwards from it with velocity V. In the 
Ishibashi-Orihara model the velocity is assumed to be 
dependent only on the electric field E, and not on the 
domain radius r(t). This makes the problem analytically 
tractable but gives rise to unphysical fitting parameters, 
such as fractional dimensionality D. The fractional D 
is not related to fractals. It is an artifact that arises 
because domain wall velocity V is actually proportional 
to l/r(t) for each domain and is not a constant at 
constant E. The volume of a domain at time t is given 
by 



C{t,t') = C D [ V(t")dt 



,1 'ii D 



(1) 



where D is the dimensionality of the growth and 
Cd is a constant which depends on the dimensional- 
ity. It is also assumed within this model that the 
nucleation is deterministic and occurs at pre-defined 
places; i.e., this is a model of inhomogeneous nucleation. 
This is an important point since some researchers still 
use homogeneous nucleation models. These are com- 
pletely inappropriate for ferroelectrics (where the nucle- 
ation is inhomo geneous, as is demonstrated by imag- 
ing experim ents fGanpule et all l200lt IShur et all l2000t 
IShurMl99fln . 

The result of the model is that the fraction of switched 
charge as a function of field and frequency may be ex- 
pressed as 



Q(£,/) = l-exp(-r D <i>(£0) 



(2) 



where $(-E) depends on the waveform used for switching. 
After some consideration and the substitution $(£') = 
E one obtains a useful relationship for the dependence 
of the field on frequency: 



For many years the standard model to de- 
scribe this process has been the Ishibashi-Orihara 



(3) 
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This r elationship has been use d to fit data fairl y well 
in TG S llrlashimoto et all I1994D. PZT and SBT (|Scotti 
199fift . More recen tly however iTsurumi et all (|200lh and 
Jung et all §002) lave found that over larger frequency 
ranges the data on several mate rials is better fi t ted b y 
the nucleation limited model of iDu and CherTI {1998). 
iTagantsev et all ()2002|) have also found that over large 
time ranges the Ishibashi-Orihara model is not a good de- 
scription of switching current data and that a nucleation 
limited model is more appropriate. It is quite possible 
of course that domain wall-limited switching (Ishibashi) 
is operative in one regime of time and field but that in 
another regime the switching is nucleation-limited. 



2. Nucleation models 

Some of the earliest detailed studies of switching in fer- 
roelectrics developed nucleation limited switching models 
where the shape of the nucleus of th e rever s ed do main 
was very important. In the work of iMer d l)l954j) and 
IWiederl l)l956j) IWiederl l)l957(l a nucleation-limited model 
was used in which when dagger-shaped nuclei were as- 
sumed, the correct dependence of the switching current 
on electric field could be derived. This approach leads to 
the concept of an activation field for nucleation (some- 
what different from the coercive field). Acti vation fields 
in thin film PZT cap acitors were measu red bv lScott et all 
(1988); very recently ljung et all l|2004|) have studied the 
effects of micro-geometry on the the activation field in 
PZT capacitors. 

The switching model of lTagantsev et all (|2002^) is a dif- 
ferent approach in which a number of non-interacting el- 
ementary switching regions are considered. These switch 
according to a broad distribution of waiting times. 



3. The scaling of coercive field with thickness 

For t he last forty years the semi - empir ical scal- 
ing law. l)Janoved 119581 iKav and DunnL Il962f) E c (d) oc 
d~ 2 f z , has been used successfully to describe the 
thickness-dependence of the coercive field in ferro electric 
films ranging from 100 microns to 20 nanometers. ij Scotd. 
|2000|) In the ultrathin PVDF films of lBune et all lll998|) a 
devi ation from this relationsh ip was seen for the thinnest 
films l)Ducharme et all l2000|) . Although they attribute 
this to a new kind of switching taking place (simultane- 
ous reversal of polari sation, as opposed to n ucleation and 
growth of domains) , iDawber et all l)2003bf) have shown, 
to the contrary, that if the effects of a finite depolar- 
ization field due to incomplete screening in the electrode 
are taken into account, then the scaling law holds over six 
decades of thickness and the coercive field does not devi- 
ate from the value predicted by the scaling law (Fig. |HJ). 
Recently iPertsev et all l)2003|) measured coercive fields in 
very thin PZT films. Although they have used a different 
model to explain their data, it can be seen that in fact 



the scaling law describes the data very well. 
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FIG. 6 The scaling of coercive field with thickness in ferr o- 
electrics; from mm to nm scale, from Daw ber~et alJ ^20031 ) 



4. Mobility of 90° domain walls 

The mobility of domain walls, especially 90-degree 
walls, depends upon their width. In this respect the 
question has been controversial, with some authors claim- 
ing very wide widths (hundreds of angstroms) and im- 
mobile walls. Some recent papers show experimentally 
that 90-degree domain walls in pero v skite f erroelectrics 
are extremely narrow llFoeth et al. I ll999D JTsai et all . 
ll992ft . l|Stemmer et all ll995|) . l)Floauet et all Il997|) . In 
PbTi03 they are 1.0 ± 0.3 nm wide. This connects the 
general question of how wide they are and wheth e r they 
are immobile. The review bv iFloauet and Valotl l)l999j) 
is quite good. They make the point that in ceramics 
these 90-degree walls are 14.0 nm wide (an order of mag- 
nitude wider than in single crystals). This could be why 
theory and experiment disagree, i.e. that something spe- 
cial in the ceramics makes them 10-15 times wider (and 
less mobile?). The latter point is demonstrated clearly 
in experiments on KNb03, together with a theoretical 
model that explains geometrical pinning in polyaxial fer- 
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roelectrics in terms of elect rostatic forces. In this respect 
the first principles study of iMever and Vanderbiltl l)2002j) 
is extremely interesting. Not only do they show that 
90 degree domain walls in PbTi03 are narrow and form 
much more easily than 180 degree domain walls, but that 
they should be much more mobile as well, the barrier for 
motion being so low they predict thermal fluctuation of 
about 12 unit cells at room temperature, which could 
perhaps explain why they appear to be wide. 

Some experimental studies using atomic force mi- 
croscopy (AFM) have attempted to answer the ques- 
tion of whether 90-degree domain walls were mobile 
or not. In certain c i rcumstances they we re immobile 
l)Ganpule e^"aZl.l2000l).llGanpule et aIl . l2000(K but in an- 
other study l)Nagaraian^7a/T ~ l2003|) the motion of 90- 
degree domain walls under an applied field was directly 
observed. It seems that in principle 90 degree domain 
walls can move, but this depends quite strongly on the 
sample conditions. 



5. Imaging of domain wall motion 

The direct imaging of ferroelectric domain walls is an 
excellent method for understanding domain wall motion 
and switching. At first this was carried out in materials 
where the domains were op tically distinct such as lead 
germanate(Sh ur et alJ . Il990|) . but more recently atomic 
force microscopy (AFM) has become a powerful tool for 
observing domain wall motion. The polarisation at a 
point can be obtained from the piezoresponse detected by 
the tip, and the tip itself can be used to apply a field to 
the ferroelectric sample and initiate switching. It is thus 
possible to begin switching events and watch their evolu- 
tion over time. AFM domain writing of ferroelectric do- 
mains can also be used to write extremely sma l l dom ain 
structures in high density arrays l|Paruch et all l200l|) or 
other d evice-like structures, such a s surface acoustic wave 
devices IjSarin Kumar et all 12004(1 . 

The backswitching studies of iGanpule et all l)200lj) 
show two very interesting effects (Fig. [TJ). The first is 
the finding that reverse domains nucleate preferentially 
at antiphase bound aries. This was stud ied in more de- 
tail subsequently bv lRoelofs et al .] (l2002t) who invoked a 
depolarisation field mediated mechanism to explain the 
result. Another explanation might be that the strain is 
relaxed at these antiphase boundaries, resulting in favor- 
able conditions for nucleation. Secondly, the influence 
of curvature on the domain wall relaxation is accounted 
for within the Kolomogorov-Avrami framework. The ve- 
locity of the domains is dependent on curvature; and as 
the relaxation proceeds, the velocity decreases and the 
domain walls become increasingly faceted. In Fig. [7| 
the white polarisation state is stable, whereas the black 
is not. The sample is poled into the black polarisation 
state and then allowed to relax back. 

A different kin d of study was undertaken by 
iTvbell et al\ l)2002[K in which they applied a voltage pulse 




FIG. 7 \ Ganmde et alJ \200A) Piezoresponse scans of a single 
cell in PbZro.2Tio.8O3, (b)-(d) illustrate the the spontaneous 
reversal of polarization within this region after wait times of 
(b) 1.01 x 10 3 , (c) 1.08 x 10 5 , (d) 1.61 x 10 s , and (d) 2.55 x 10 5 
s. Faceting can be seen in (c),(d),and (e). (f) Transformation- 
time curve for the data in (b)-(e) 

to switch a region of the ferroelectric using an AFM tip 
and watched how the reversed domain grew as a function 
of pulsewidth and amplitude. They were able to show 
that the process was well described by a creep mech- 
anism, thought to arise due to random pinning of do- 
main walls in a disordered system (Fig. (SJ. Though 
the exact origin of the disorder was not clear, it is sug- 
gested that it is connected to oxygen vacancies, which 
also play a role in pinning the domai n walls during a 
fatigue process due to their or dering ijPark and Chadil 
1998h l|Scott and Dawbe"rll2000j) . 
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Pulse width (s) i/Edo V'cm) 

FIG. 8 \Tvbell et all \200A) (a) Domain size increases loga- 
rithmically with pulse widths longer than 20 /xs and saturates 
for shorter times as indicated by the shaded area, (b) Domain 
wall speed as function of the inverse applied electric field for 
290, 370, and 810 A thick samples. The data fit well the 
characteristic velocity field relationship of a creep process. 

An unsolved puzzle is the direct observa- 
tion via AFM of domain walls penetrating grain 
boundaries ljGruverman et all Il997fl . This is contrary 
to some expectations and always occurs at non-normal 
incidence, i.e. at a small angle to the grain boundary. 

One of the interesting things to come out of the work in 
lead germanate (where ferroelectric domains are optica lly 
distinct due to electrogyration) bv IShur et all 1|1990|) is 
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that at high applied electric fields (15kV cm -1 ) tiny do- 
mains are nucleated in front of the moving domain wall 
(Fig. EJ. A very similar effect is seen in ferromagnets as 
observed by iRandoshkinl (^95) in a single crystal iron 
garnet film. 




FIG. 9 \Shur et alJ 1 199(1) Nucleation of nanodomains in front 
of domain wall in lead germanate at high electric field. Black 
and white represent the two directions of polarisation. 



However in ferroma gnets the effect is modelled by a 
spin wave mechanism. ijKhodenkovt Il975(l This mecha- 
nism is based on the gyrot ropic model of domain wall mo- 
tion in uniaxial materials! Walkcil ll963|) . When a strong 
magnetic driving field (exceeding the Walker Threshold) 
acts upon a domain wall, the magnetization vectors in 
the domain wall begin to precess with a frequency jH, 
where 7 is the effective gyromagnetic ratio. By relat- 
ing the precession frequency in the domain wall with the 
spin wave frequency in the domain, good predictions can 
be made for the threshold fields at which the effect oc- 
curs. We note that the domains nucleated in front of the 
wall may be considered as vortex-like skyrmions. The 
similarity between these effects is thus quite surprising 
and suggests that perhaps there is more in common be- 
tween ferroelectric domain wall motion and ferromagent 
domain wall motion than is usually c onsidered. How- 
ever, whereas fDemocritov et al. 1 l)l988|) have shown that 
magnetic domain walls can be driven supersonically (re- 
sulting in a phase-matched Cerenkov-like bow wave of 
acoustic phonon emission), there is no direct evidence 
of supersonic ferroelectric domains. Processes such as 
the nano-domain nucleation described above seem to oc- 
cur instead when the phase velocity of the domain wall 
motion approaches the speed of sound. Of course the 
macroscopic electrical response to switching can arrive 
at a time t < v/d where v is the sound velocity and d, 
the film thickness, simply from domain nucleation within 
the interior of the film between cathode and anode. 



B. Electrical Characterization 

1. Standard Measurement techniques 

Several kinds of electrical measurements are made on 
ferroelectric capacitors. We briefly introduce them here 



before proceeding to the following chapters where we dis- 
cuss in detail the experimental results obtained by using 
these techniques. 

a. Hysteresis One of the key measurements is naturally 
the measurement of the ferroelectric hysteresis loop. 
There are two measurement schemes commonly used. 
Traditionally a cap acitance bridge as fi rst de scribed by 
Sawyer and Tower l)Sawver and Toweit ^30') was used 
(Fig. 111))) . Although this is no longer the standard 
way of measuring hysteresis the circuit is still useful 
(and very simple and cheap) and we have made sev- 
eral units which are now in use in the teaching labs 
in Cambridge for a demonstration in which students 
are able t o make and test their own ferroelectric KNO3 
capacitor l|Dawber et all l2003a|) . 




IS- It" -IT 



(a) (b) 

FIG. 10 \Sawver and Towe~k \l9S(\) fa)The original Sawyer- 
Tower circuit, (b) Hysteresis in Rochelle salt measured using 
this circuit by Sawyer and Tower at various temperatures 

This method is not very suitable in practice for 
many reasons, for example the need to compensate 
for dielectric loss and the fact that the film is be- 
ing continuously cycled. Most testing of ferroelec- 
tric capacitors is now carried out using commercial 
apparatus from one of two companies, Radiant Tech- 
nologies (http://www.ferrodevices.com/) and AixAcct 
(http://www.aixacct.com/). Both companies' testers 
can carry out a number of tests and measurements, and 
both machines use charge or current integration tech- 
niques for measuring hysteresis. Both machines also offer 
automated measurement of characteristics such as fatigue 
and retention. 

In measuring P(E) hysteresis loops several kinds of ar- 
tifacts can arise. Some of these are entirely instrumen- 
tal, and some arise from the effects of conductive (leaky) 
specimens. 

Hysteresis circuits do not measure polarization P di- 
rectly. Rather, they measure switched charge Q. For an 
ideal ferroelectric insulator 

Q = 2P r A (4) 

where P r is the remanent polarization and A is elec- 
trode area is a parallel plate capacitor. For a somewhat 
conductive sample 
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Q = 2P r A + <jE a t (5) 

where a is the electrical conductivity; E a , is the applied 
field and t, the measuring time. Thus Q in a pulsed 
measuring system depends on the pulse width. 

The four basic types of apparent hysteresis curves that 
are artifacts are shown in Figllll 




(Breakdown) 



FIG. 11 Common hysteresis artifacts: (a) Dead short, (b) 
linear lossy dielectric, (c) saturated amplifier, (d) non-linear 
lossy dielectric. 

Fi dlll a is a dead short in a Sawyer Tower circuit or 
modern variant and is discussed in the instruction docu- 
mentation for both the AixAcct 1 and Radiant 2 testers. 

Fig ^] b shows a linear lossy dielectric. The points 
where the loop crosses V a = are often misinterpreted 
as P r values. Actually this curve is a kind of Lissajous 
figure. It can be rotated out of the page to yield a 
straight line (linear dielectric response). Such a rota- 
tion can be done electrically and give a "compensated" 
curve. Here "compensation" means to compensate the 
phase shift caused by dielectric loss. 

Fig 1111 c is more subtle. Here are two seemingly per- 
fect square hysteresis loops, obtained on the same or 
nominally equivalent specimens at different maximum 
fields. The smaller loop was run at an applied voltage 
of V a = 10V, and yields P r = QiLCcmT 2 and the larger 
at V a = 50V and yields P r = lOO^Ccm -2 . Note that 
both curves are fully saturated (flat tops). This is im- 
possible, if the dipoles of the ferroelectric are saturated 
at P r — 6fiCcm~ 2 then there are no additional dipoles 
to produce P r = 100/iCcm~ 2 in the larger loop at high 
voltage. What actually occurs in the illustration is sat- 
uration of the amplifier in the measuring system, not 
saturation of the polarization in the ferroelectric. The 



1 TF Analyzer 2000 FE-module instruction manual 

2 http:/ /www. ferroelectrictesters. com/html/specs. html#:tut 



figure is taken from I.Taffe. Cook and Jaffe 1 119711) where 
this effect is discussed (p. 39). It will be a serious prob- 
lem if conductivity is large in Eq. 03 "Large" in this sense 
is a > lO^iflcm)- 1 and "small" is <r < 10 _7 (f2cm) -1 . 
This is probably the source of P r > \hO[iCcmT 2 reports 
in BiFeC>3 where a can exceed 10 -4 (Q.cm)^ 1 . 

Finally Fig 1111 is a nonlinear lossy dielectric. If it is 
phase compensated it still resembles real hysteresis. One 
can verify whether it is real or an artifact only by varying 
the measuring frequency. Artifacts due to dielectric loss 
are apt to be high ly frequency dependent . Figs ^2 b and 
d are discussed in iLines and Glass] lll977l) (p.104). 

No data resembling Figs llll a-d should be published as 
ferroelectric hysteresis. 



b. Current measurements Another measurement of im- 
portance which is carried out in an automated way by 
these machines is the measurement of the leakage cur- 
rent. This is normally discussed in terms of a current- 
voltage (I-V) curve, where the current is measured at a 
specified voltage. It is important however that sufficient 
time is allowed for each measurement step that the cur- 
rent is in fact true steady state leakage current and not 
relaxation current, and for this reason current-time (I-t) 
meas urements can also be important ( Dietz and Waseit 
119951) . Relaxation times in ferroelectric oxides such as 
barium titanate are typically 1000 s at room tempera- 
ture. 



c. Dielectric permittivity An impedance analyser mea- 
sures the real and the imaginary parts of the impedance 
by use of a small-amplitude AC signal which is applied 
to the sample. The actual measurement is then made 
by balancing the impedance of the sample with a set of 
references inside the impedance analyser. From this the 
capacitance and loss can be calculated (all this is done 
automatically by the machine). It is possible at the same 
time to apply a DC bias to the sample, so the signal is 
now a small AC ripple superimposed on a DC voltage. 
Ferroelectric samples display a characteristic "butterfly 
loop" in their capacitance voltage relationships, because 
the capacitance is different for increasing and decreas- 
ing voltage. The measurement is not exactly equivalent 
to the hysteresis measurement. In a capacitance voltage 
measurement a static bias is applied and the capacitance 
measured at that bias, whereas in a hysteresis measure- 
ment the voltage is being varied in a continuous fashion. 
Therefore it is not strictly true that C(V) = (d/dV)P(V) 
[area taken as unity] as claimed in some texts, since the 
frequency is not the same in measurements of C and 
P. Impedance spectroscopy (where the frequency of the 
AC signal is varied) is also a powerful tool for analy- 
sis of films, especially as it can give information on the 
timescales at which processes operate. Interpretation of 
these results must be undertaken carefully, as artifacts 
can arise in many circumstances and even when this is 
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not the case many elements of the system (e.g., elec- 
trodes, grain boundaries, leads etc.) can contribute to 
the impedance in complicated ways. 



2. Interpretation of dielectric permittivity data 

a. Depletion charge vs intrinsic response Before looking 
for ferroelectric contributions to a system's electrical 
properties one should make sure there are not contri- 
butions due to the properties of the system unrelated 
to ferroelectricity. Although much is sometimes made 
of the dependence of capacitance on voltage, it is worth 
noting that metal-semiconductor-metal systems have a 
characteristic capacitance voltage which arises from the 
response of depletion layers to applied voltage (Fig. H2)l . 
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FIG. 12 Current voltage and capacitance voltage re- 
lation ship of Pt-Si-Pt punch-through diode \Sze et all 
\197A) , characteristics very similar to those obtained in metal- 
ferroelectric-metal systems 



The change in with field due to its nonlinea ri ty ha s also 
been calculated b y both lOutzourhit et all l|l995|) and 
iDietz et all l)l997j) . The real problem is that both pic- 
tures are feasible. One should not neglect the fact that 
the materials have semiconductor aspects; but at the 
same time it is not unreasonable to expect that the known 
nonlinearity of the dielectric response in these materials 
should be expressed in the capacitance voltage character- 
istic. Probably the best approach is to avoid making any 
conclusions on the basis of these kinds of measurements 
alone, as it is quite possible that the relative sizes of the 
contributions will vary greatly from sample to sample, 
or even in the same sample under different experimental 
conditions. 
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FIG. 13 ( Capa citance v applied bias for BST thin films. The 
data are due to \Basceri et alJ \199'\) 



Essentially the problem boils down to the fact that 
there are two possible sources of the dependence of ca- 
pacitance on applied field, either changes in depletion 
width or changes in the dielectric constant of the mate- 
rial, i.e. 



C(E) e(E) 



A 



d(E) 



(6) 



Several groups have assumed, that all the change in 
the capacitance with field C(E) comes from change in 
depletion width d(E) and that (E) is changing negligi- 
bly. The first to suggest this was Joe Evans (IEEE- 
ISAF Meeting, Urbana, 111. 1990), who found d = 20 
nm in PZT. Later Sandia cl aimed that there w as no de- 
pletion (d=0 or d=innnitv) <|Mirler et all . Il990(l . Several 



b. Domain wall contributions Below the coercive 
field there are also contributions to the permittivity 
from dom ain walls, as first p ointed out by Fouskova 
in 1965 l)Fouskoval ll965j) . l)Fouskova and Janousekl 
1965). In PZT the contributions of domain 
wall pinning to the dielectric permittivity has 
been studied in d etail by Damianovic and 
Taylo r l|Damianovj cl Il997t) JTavlor and Damianovid 
ll997l) . lfravlor and Damianovicl Il998|) . who showed that 
the sub-coercive field contributions of the permittivity 
were described by a Raleigh law with both reversible 
and irreversible components, the irreversible component 
being due to domain wall pinning. 



(iBrennanl .19921). llMihara et all .1992). (Hwane et al. 


1998), ijHwand. 


Il998|) .((Scott et all. 


11992(1. ((Saver et al. 



In contrast to the approach of explaini ng these char- 
acteri stics using semiconductor models iBasceri et all 
(1199711 account for their results on the basis of a Landau- 
Ginzburg style expansion of the polarisation (Fig. I13fl . 



c. Dielectric measurements of phase transitions One of the 
most common approaches to measuring the transition 
temperature of a ferroelectric material is naturally to 
measure the dielectric constant and loss. However, in 
thin films there are significant complications. In bulk 
the maximum in the dielectric constant is fairly well cor- 
related with the transition temperature, but this does not 
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al ways seem to be the case in thin films. As pointed out 
by IVendik and Zubkol l)2000|) . a series capacitor model 
is required to extract the true transition temperature, 
which in the case of BST has been show n to be indepen- 
dent of thickness ljLookman et al\ . |2004|) . in contrast to 
the temperature at which the permittivity maximum oc- 
curs, which can depend quite strongly on thickness (Fig. 
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FIG. 14 lLookman et alJ 1200A) Comparison between the ap- 
parent Curie temperature in BST taken from Curie-Weiss 
plots of raw data (empty circles) and intrinsic data after cor- 
rection for interfacial capacitance (filled squares) had been 
performed. The intrinsic Curie temperature appears to be in- 
dependent of film thickness 



3. Schottky barrier formation at metal-ferroelectric junctions 

In general, since ferroelectric materials are good insu- 
lator the majority of carriers are injected from the elec- 
trode. When a metal is attached to a ferroelectric ma- 
terial, a potential barrier is formed if the metal work 
function is greater than the electron affinity of the ferro- 
electric. This barrier must be overcome if charge carriers 
are to enter the ferroelectric, On the other hand if the 
electron affinity is greater than the work function, then 
an ohmic contact is formed. For the usual applications 
of ferroelectrics (capacitors) it is desirable to have the 
largest barrier possible. If a metal is brought into contact 
with an intrinsic pure ferroelectric and surface states do 
not arise (i.e., the classic metal- insulator junction), then 
the barrier height is simply 



X 



(7) 



In this case the fermi level of the metal becomes the 
fermi level of the system, as there is no charge within the 
insulator with which to change it. On the other hand, 
if there are dopants or surface states, then there can be 
a transfer of charge between the metal and ferroelectric, 
which allows the ferroelectric to bring the system fermi 
level towards its own fermi level. 



Although single crystals of undoped ferroelectric ti- 
tanates tend to be slightly p-type, simply because there 
are greater abundances of impurities with lower valences 
than those of the ions for which they substitute (Na+ 
for Pb +2; Fe+3 f or Ti+4) llSmvtbi ll984h . llChan et all 
ll98lD . llChan et 

aD,ES76) m reality most ferroelectric ca- 
pacitors are fine grained polycrystalline ceramics and are 
almost always oxygen deficient. Typically the regions 
of the capacitor near to interface are more oxygen de- 
ficient than the bulk. Oxygen vacancies act as donor 
ions, and this means there can be a transition from n- 
type behaviour at the interface to p-type behaviour in 
the centr e of the film as is ev i dent i n the kelvin probe 
study of iNowotnv and Rekas I l)l994|) , who found that 
in bulk BaTiOa with Pt electrodes a change in work- 
function from 2.5 ± 0.3eV for surfaces and 4.4 ± 0.4eV 
in the bulk of the material. The nature of the material 
near the surface is important since it determines whether 
a b locking or ohmic contact i s formed. It has been shown 
bv lDawber and Scotfi J2002) that the defect concentra- 
tion profile as measured by capacitance voltage technique 
may be explained by a model of combined bulk and grain 
boundary diffusion of oxygen vacancies during the high 
temperature processing of a film. 

Regardless of the p-type or n-type nature of the mate- 
rial , in most oxide ferroelectrics on elemental metal elec- 
trodes the barrier height for electr ons is significantly less 
than the barrier height for holes ljRobertson and ChenL 
119991) . and so the dominant injected charge carriers are 
electrons. As the injected carriers dominate the conduc- 
tion, leakage currents in ferroelectrics are electron cur- 
rents and not hole currents, contra ry to the suggestion of 
Stoli chnov and Taeantsevl l)l998bj) . 

The first p icture of barrier form ation i n sem iconduc- 
tors is due to lSchottkvl <jl93S|) and lMotd <jl93Sft . In this 
picture the conduction band and valence band bend such 
that the vacuum levels at the interface are the same and 
the fermi level is continuous through the interface, but 
deep within the bulk of the semiconductor retains its orig- 
inal value relative to the vacuum level. This is achieved 
by the formation of a depletion layer which shifts the 
position of the fermi level by altering the number of elec- 
trons within the interface. 

Motivated by the experimental observation that many 
Schottky barrier heights seemed to be fairly indepen- 
dent of the metal used for the electrode iBardeenl lll947ft 
proposed a different model of metal-semiconductor junc- 
tions. In this picture the fermi level of the semiconductor 
is "pinned" by surface states to the original charg e neu- 
trality level. These states, as first suggested by iHeind 
(1965), are not typically real surface states but rather 
states induced in the band-gap of the semiconductor by 
the metal. 

Most junctions lie somewhere between the Schottky 
and Bardeen limits. The metal induced gap states 
(MIGS) can accommodate some but not all of the differ- 
ence in the fermi level between the metal and the semi- 
conductor, and so band bending still occurs to some ex- 
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tent. The factor S — is used to define this, with 
S = 1 being the Schottky limit and S = represent- 
ing the Bardeen limit. The value of S is determined by 
the nature of the semiconductor; originally experimen- 
tal trends linking this to the covalency or ionicity of the 
bonding in the material were observed (with covalent ma- 
teria ls developing many more MIGS. than ionic materi- 
als) ijKurtin et aZ.l . ll96£y . However better correlation was 
found between the effective band g aps (dependent on the 
electronic dielectric constant f m ) l|Schluteii[T97^ . with 
llMonchLlT98l 



(--l) = 0.1( £oo -l) 2 . 



(8) 



Although SrTiC>3 was invoked as one of th e materi- 
als th at violated the electronegativity rule by ISchluterl 
||1978|) . it is omitted from the plot against (e^ — 1). 
The experimental value for S in SrTiC>3 can be measured 
from Dietz's data as approximately 0.5. This does not 
agree well with what one would expect from Monch's 
empirical relation, which g ives S = 0.28 (as used by 
iRobertson and Chenl (^999)). Note that the use of the 
ionic trap-free value S=l for BST gives a qualitative er- 
ror: It predicts that BST on Al should be ohmic, whereas 
in actuality it is a blocking junction; an S-value of ap- 
proximately 0.3 predicts a 0.4 e V Schottky b arrier height, 
in agreement with experiment. (jScottL l2000f ) 

We can extract the penetration depth for Pt states into 
BaTiO .s from the first principles calculation of lRao et all 
i|l997j) by fitting the density of states (DOS) of platinum 
states in the oxygen layers to an exponential relationship 
t o extract the characteris tic length as 1.68 angstroms. 

ICowlev and S^el <ll965l) derived an expression for the 
barrier height for junctions between the two extremes. 
In this approach the screening charges in the electrode 
and the surface states are treated as delta functions of 
charge separated by an effective thickness <5 e / / . This ef- 
fective thickness takes into account both the Thomas- 
Fermi screening length in the metal and the penetra- 
tion length of the MIGS, and is essentially an air-gap 
approach. 

The expression for the barrier height is 



S 2 C 



h = S{cp m - X ) + (1 - S)(E g - <M + C 
Si[C( ( / )m -x) + (l-S)(E g -4, Q )^ 



(9) 



^(E g -E f + kT) 



(10) 



In the above S 
s « lOeo an d No 



,C 



2qe s N D S 2 eff . When 
3 C is of the order of 



e s PS We and 1\ D < 10 18 cm 
0.01 eV and it i s reas onable to discard the term £ as 
ICowlev and Szd l)l965[) did. Neglecting this term, a s 
has been pointed out bv lRhoderick and Williams] l)l988[) . 
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FIG. 15 Energy band diagram of a metal n-type semiconduc- 
tor contact after Cowley and Sze. 



is equivalent to neglecting the charge in the depletion 
width. In the systems under consideration here this term 
should not be neglected as it can be quite large. To 
demonstrate the effect on the barrier height we calculate 
the barrier height for a Pt-SrTi03 barrier over a wide 
range of vacancy concentrations (Fig. 116(1 . 
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FIG. 16 Schottky barrier height of Pt-SrTiOz as a function 
of oxygen vacancy concentration, note that this may explain 
the variation of experimental values from ca. 0.7 eV to 1.0 
eV 

It can be seen that the effect of vacancies on bar- 
rier height becomes important for typical concentra- 
tions of vacan c ies en countered in ferroelectric thin films. 
iDawber et al\ l|200l|) have addressed this issue and also 
the effect of introduced dopants on barrier heights. De- 
spite their omission of the term discussed above the work 
of Robertson and Chen is valuable because of their calcu- 
lation of the charge neutrality levels for several ferroelec- 
tric materials, an essential parameter for the calculation 
of metal- ferroelectric barrier heights. 

In a ferroelectric thin film this distribution of charges 
at the interface manifests itself in more ways than sim- 
ply in the determination of the Schottky barrier height. 
Electric displacement in the system is screened over the 
entire charge distribution. 

In measuring the small -signal capacitance against 
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thickness there is always a non-zero intercept, which has 
been typically associated with a "dead layer" at the metal 
film interface. However, in most cases this interfacial 
capacitance can be understood by recognizing that a fi- 
nite potential exists across the charge at the interface. 
In the simplest approximation one neglects any charge 
in the ferroelectric and uses a Thomas-Fermi screening 
model for the metal. This was initially considered by 
iKu and UllmarJ ll 19641) and f i rst ap plied to high k di- 
electrics bv lBlack and Welserl ||1999(1 . In their work they 
use a large value for the dielectric constant of the oxide 
metal, considering it as the dielectric response of the ions 
stripped of their electrons. This may seem quite reason- 
able but is not however appropriate. In general we think 
of metals not being able to sustain fields, and in the bulk 
they certainly cannot, but the problem of the penetra- 
tion of electric fields into metals is actually well known 
in a different context, that of the microwave skin depth. 
It is very instructive to go through the derivation as an 
AC current problem and then find the DC limit which 
will typically apply for our cases of interest. 

We describe the metal in this problem using the Drude 
Free Electron Theory: 



1 + IU1T 



ill) 



There are three key equations to describe the charge 
distribution in the metal. 

Poisson's equation for free charges: 



The continuity equation: 



1 8E(z) 
47r dz 



iujp{z) 



dj{z) 
dz 



The Einstein transport equation: 
These are combined to give: 

8 2 p 47TCT ILU 



(12) 



(13) 
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This tells us that if at a boundary of the metal there 
exists a charge it must decay with the metal exponen- 
tially with characteristic screening length A: 



(16) 



In the DC limit (which applies for most frequencies of 
our interest) this length is the Thomas-Fermi screening 
length; 



,47T0£ 

1 D ' 



-1/2 



So it becomes clear that the screening charge in the 
metal may be modelled by substituting a sheet of charge 
displaced from the interface by the Thomas- Fermi screen- 
ing length, but that in calculating the dielectric thickness 
of this region the effective dielectric constant that must 
be used is 1, consistent with the derivation of the screen- 
ing length. Had we used a form of the Poisson equation 
that had a non-unity dielectric constant, i.e., 



p(z) 



e dE{z) 



An dz 

then our screening length would be 



(18) 
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which is not the Thomas-Fermi screening length. Thus 
the use of a non-unity dielectric constant for the metal is 
not compatible with the use of the Thomas- Fermi screen- 
ing length. 

Measurements on both sol-gel and CVD lead zirconate- 
titanate (PZT) films down to ca. 60 nm thickness show 
that reciprocal capacitance 1/C(d) versus thickness d ex- 
trapolates to finite values at d=0, demonstrating an in- 
terfacial capacitance. However whereas the value for the 
sol-gel films is consistent with the Thomas Fermi screen- 
ing approach (.05nm), the value of interfacial thickness 
(0.005 nm) for the CVD films is only 10% of the in- 
terfacial capacitance that would arise from the known 
Fermi- Thomas screening leng th of 0.05 nm in the Pt elec- 
trodes l)Dawber et all 12003 b|) . That is, if this result were 
interpreted in terms of a "dead layer", the dead layer 
would have negative width. This result may arise from a 
compensating "double-layer" of space charge inside the 
semiconducting PZT dielectric; t he Armstrong- Horrocks 
((Armstrong and HorrocksT Il997[) semiconductor formal- 
ism form of the earlier Helmholtz and Gouy-Chapman 
polar-liquid models of the double layer can be used. Such 
a double layer is unnecess ary in PVDF be cause that ma- 
terial is highly insulating l|Moreira 1 12002'). This explains 
quantitatively the difference (x8) of interfacial capaci- 
tance in sol-gel PZT films compared with CVD PZT films 
of the same thickness. The magnitude of the electroki- 
netic potential (or zeta-potential) C = f^" that develops 
from the Helmholtz layer can be estimated without ad- 
justable parameters from the oxygen vacancy gradient 
data of Dey for a typical oxide perovskite, SrTi03; us- 
ing Dey's surface charge density a of 2.8 x 10 18 e / m 2 , 
a Gouy screening length in the dielectric d' = 20 nm, 
and a dielectric constant of e = 1300 yields C = 0.78 eV; 
since this is comparable to the Schottky barrier height, it 
implies that much of the screening is provided internally 
by mobile oxygen vacancies. [Here cr(r, fi) is a function 
of time r and mobility fi for a bimodal (a.c.) switching 
process.] 
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4. Conduction mechanisms 

In general conduction is undesirable in memory devices 
based on Capacitors, and so the understanding and min- 
imization of conduction has been a very active area of 
research over the years. Many mechanisms have been 
proposed for the conduction in ferroelectric thin films. 



a. Schottky injection Perhaps the most commonly ob- 
served currents in ferroelectrics are due to thermionic in- 
jection of electrons from the metal into the ferroelectric. 
The current-voltage characteristic is determined by the 
image force lowering of the barrier height when a poten- 
tial is applied. A few points should be made about Schot- 
tky injection in ferroelectric thin films. The first is about 
the dielectric constant appropriate for use. In ferro- 
electrics the size of the calculated barrier height lowering 
depends greatly on which dielectric constant, the static or 
the electronic, is used. The correct dielectric constant is 
the electronic one (~ 5. 5), as discussed by Scott and used 
by Dietz, and by Zafar. iDietz an cLffiasgj fcLMjj) u s ed the 
more general injection law of iMurnhv and Good! l|l956h 
to describe charge injection in SrTiOa films. They found 
that for lower fields the Schottky expression was valid, 
but at higher fields numerical calculations using the gen- 
eral injection law were required. They did not however 
find that Fowler-Nordheim Tunneling was a good descrip- 
tion of any of the experimental data . 

It has been shown by Za far et~all l)l998[) that in fact 
the correct form of the Schottky equation that should 
be used for ferroelectric thin films is the diffusion lim- 
ited equation of S i mmon s. Furthermore, very recently 
iDawber and Scottl l)2004f) have shown that when one con- 
siders the ferroelectric capacitor as a metal-insulator- 
metal system with diffusion-limited current (as opposed 
to a single metal- insulator junction), the leakage cur- 
rent is explained well; in addition, a number of unusual 
effects,, su ch as the negative diff erential resistivity ob- 
served by IWatanabe et all Jl998h and the PTCR e ffect 
observed bv lHwang et all l)l997ftlHwang et all l)l998f) are 
accounted for. 



b. Poole-Frenkel One of the standard ways of identifying 
a Schottky regime is to plot log(^) against V'. In this 
case the plot will be linear if the current injection mech- 
anism is Schottky injection. Confusion can arise because 
carriers can also be generated from internal traps by the 
Poole-Frenkel effect, which on the basis of this plot is 
indistinguishable from Schottky injection. However, if 
the I-V characteristic is asymmetric with respect to pos- 
itive and negative voltages (as is usually the case) then 
the injection process is most probably Schottky injection. 
There are however some papers that show symmetrical 
I-V curves and correctly explain their data on the basis 
of a Poole-Frenkel conduction mechanism l|Chen et all . 
1998). 



c. Fowler-Nordheim Tunneling Many researchers have 
discussed the possibility of tunnelling currents in ferro- 
electric thin film capacitors. For the most part they are 
not discussing direct tunneling through the film, which 
would be impossible for typical film thicknesses, but in- 
stead tunneling through the potential barrier at the elec- 
trode. The chief experiment al evidence that it might in- 
deed be possible is due to IStolichnov. Taeantsev et all 
(1998) who have seen currents that they claim to be 
entirely tunneling currents in PZT films 450 nm thick 
at temperatures between 100-140 K. It should be noted 
however that they only observed tunneling currents 
above 2.2MV/cm, below which they were unable to ob- 
tain data. The narrowness of the range of fields for which 
they have collected data is a cause for concern, since the 
data displayed in their paper go from 2.2-2.8 MV/cm. 
We conducted leakage current measurements on a 70 nm 
BST thin film at 70 K and found that the leakage cur- 
rent, while of much lower magnitude, was still well de- 
scribed by a Schottky injection relationship, although if 
one fitted this data to a similarly narrow field region it 
did appear to satisfy the Fowler-Nordheim relationship 
well (Fig. El • 
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FIG. 17 Leakage current data from Au-BST-SrRuOs film at 
room temperature and at T=70K 

The effectiv e masses for tunneling obtaine d in 
the studies of Stolichnov, Taeantse v et al. | II1999D and 
iBaniecki et all l|200lfl also seem to be at odds with 
the normal effective masses considered for these mate- 
rials. Whereas they use effective masses of 1.0, the ef- 
fective masses in perovskite oxides seem to be some- 
what larger, fr om 5-7 m P f o r bari um titanate and stron- 
tium titanate. l)Scott et all . [2003). Although the tun- 
neling mass and the effective (band) mass need not be 
the same in general, if the tunneling is through thick - 
nesses of > 2 nm, they ar e near lyso. jSchnurml ill 967ft: 
IConlev and Mahianl l)l967(l . also find that the tunnelling 
mass due to light holes in GaAs fits the band mass very 
well. ]." 

d. Space Charge Limited Currents The characteristic 
quadratic relationship between current and voltage that 
is the hallmark of space charge limited currents are often 



16 



seen in ferroelectrics. Sometimes it is observed that space 
charge limited currents are seen when a sample is biased 
in one direction, whereas for the opposite bias Schottky 
injection dominates. 



e. Ultra-Thin Films - Dire c t Tun neling Very recently 
iRodriguez Contreras et a/1 ()2003[) have succeeded in 
producing metal-PZT-metal junctions sufficiently thin 
(6nm) that it appears that direct tunneling or phonon 
assisted tunneling (in contrast to Fowler Nordhcim tun- 
neling)through the film may occur, though this result 
requires more thorough investigation, since the authors 
note the barrier heights extracted from their data using 
a direct tunelling model are much smaller than expected. 
The principal result of this paper is resistive switching, 
which may be of considerable interest in device appli- 
cations, but also requires more thorough investigation. 
This very interesting experimental study raises impor- 
tant questions about the way that metal wave functions 
penetrating from the electrode and ferroelectric polarisa- 
tion interact with each other in the thinnest ferroelectric 
junctions. 



f. Grain boundaries Grain boundaries are often consid- 
ered to be important in leakage current because of the 
idea that they will provide conduction pathways through 
the film. 

Gruverman's results suggest that this is not the case in 
SBT. In his experiment an AFM tip is rastered across the 
surface of a polycrystalline ferroelectric film. The imaged 
pattern records the leakage current at each point: white 
areas are high-current spots; dark areas, low current. If 
the leakage were predominantly along grain boundaries, 
we should see dark polyhedral grains surrounded by white 
grain boundaries, which become brighter with increasing 
applied voltage. In fact, the opposite situation obtains: 
This indicates that the grains have relatively low resis- 
tivity, with high-resistivity grain boundaries. The second 
surprise is that the grain conduction comes in a discrete 
step; an individual grain suddenly "turns on" (like a light 
switch). Smaller grains generally conduct at lower volt- 
ages (in accord with Maier's theory of space charge effects 
being larger in small grains w ith higher surface-volume 
ratios l|Lubomirskv et all . 12002^ ) . 



C. Device Failure 

1. Electrical Breakdown 

The process of electrical shorting in ferroelectric PZT 
was first shown by Plumlee (1967) to arise from dendrite- 
like conduction pathways through the material, initiated 
at the anodes and/or cathodes. These were manifest as 
visibly dark filamentary paths in an otherwise light ma- 
terial when viewed through an optical microscope. They 



have been thought to arise as "virtual cathodes" via the 
growth of ordered chains of oxygen-deficient material. 
This mechanism was modeled in detail b y Duike r et al 
llDuiker lll990|) . l)Duiker and Beale Lll990|) . l)Duiker et all 
1199ft) . 

To establish microscopic mechanisms for beakdown in 
ferroelectric oxide films one must show that the depen- 
dences of breakdown field upon film thickness d, 
ramp rate, temperature, doping, and electrodes are sat- 
isfied. The dependence for PZT upon film thickness is 
most compatible wit h a low power-law dependence or 
possibly logarithmic ijScott et al I l2003|) . The physical 
models compatible with this include avalanche (logarith- 
mic) , collision ionization fro m electrons injected via fiel d 
emission from the cathode l|Forlani and Minna ial ll964(K 
which gives 

E B = Ad- W (20) 

with j < w < - L or the linked defect model of 
iGerson and Marshalll (|1959T) . which has d = 0.3. The 
dependence on electrode material arises from the elec- 
trode work function and the ferroelectric electron affinity 
through the r esulta nt Schott ky bar r ier he ight. Following 
IVon Hipoel 1 1119351) we have IScottl ll200fl) . p62.) 

eE B X = h($ M - (21) 

where $m and <&fe are the work functions of the metal 
and of the semiconducting ferroelectric,; is electron mean 
free path; and h is a constant of order unity. 

Even in films for which there is considerable Poole- 
Frenkel limitation of current (a bulk effect), the Schot- 
tky barriers at the electrode interfaces will still dominate 
breakdown behavior. 

In general electrical breakdown in ferroelectric oxides 
is a hybrid mechanism (like spark discharge in air) in 
which the initial phase is electrical but the final stage 
is simple thermal run-away. This makes the dependence 
upon temperature complicated. 

There are at least three different contributions to the 
temperature dependence. The first is the thermal prob- 
ability of finding a hopping path through the material. 
Following Gerson and Marshall and assuming a random 
isotropic distribution of traps, Scott l)Scottlll995|) showed 
that 



E B =G-^-\ogA (22) 

which gives both the dependence on temperature T 
and electrode area A in agreement with practically all 
experiments on PZT, BST, and SBT. 

In agreement with this model the further assumption of 
exponential conduction (non-oh mic) estimat ed to occur 
for applied field E > 30 MV/m (IScottl l2000l) 
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a(T) = a Q cxp — - (23) 

in these materials yields the correct dependence of 
breakdown time ts upon field 

\ogt B =ci -c 2 E B (24) 

as well as the experimentally observed dependence of 
Eg on rise time t c of the applied pulse: 

E B =c 3 t~ i (25) 

Using the same assumption of exponential conduction, 
which is valid for 



12911 . Even if the ferroelectrics were single crystals, with 
10 20 cm~ 3 oxygen vacancies near the surface, any In- 
dependence from Eq. [2H would be unmeasurably small; 
and for the actual fine-grained ceramics (40 nm grain di- 
ameters), the mean free path is ca. 1 nm and limited by 
gr ain boundaries (T-indep endent). Thus the conclusion 
of IStolichnov et all ()200l|) regarding avalanche is quali- 
tatively and quantitatively wrong in ferroelectric oxides. 



2. Fatigue 

Polarisation fatigue, which is the process whereby the 
switchable ferroelectric polarisation is reduced by repet- 
itive electrical cycling is one of the most serious device 
failure mechanisms in ferroelectric thin films. It is most 
commonly a problem when Pt electrodes, desirable be- 
cause of their high work-functions, are used. 



aeE <C KbT 



(26) 



where a is the lattice nearest neighbor oxygen-site hop- 
ping distance (approximat ely a lattice constant) and e, 
the electron charge, Scott llScottl. i000) shows that the 
general breakdown field expression 
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(27) 



in the impulse approximation (in which the second 
term in the above equation is neglected) yields 
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FIG. 18 Change in the polarisation hysteresis loop with fa- 
tiaue \Scott and Poulianil tl98&) 
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which suffices to estimate the numerical value of break- 
down field for most ferroelectric perovskite oxide films; 
values approximating 800 MV/m are predicted and mea- 
sured. 

A controversy has arisen regarding the temperature 
dependence of En(T) and the possibility of avalanche 
l|Stolichnov et all l200lj) In low carrier concentration sin- 
gle crystals, especially Si, avalanche mechanisms give a 
temperature dependence that is controlled by the mean 
free path of the injected carriers. This is physically be- 
cause at higher temperatures the mean free path A de- 
creases due to phonon scattering and thus one must apply 
a higher field Eb to achieve avalanche conditions. 



A = A tanh(^) 



(29) 



However, this effect is extremely small even for low 
carrier concentrations (10% change in Eb between 300K 
and 500K for n = 10 16 cm -3 ) and negligible for higher 
concentrations. The change in Eb in BST between 600K 
and 200K is > 500% and arises from Eq. [Z2\ not Eq. 



Importantly fatigue occurs through the pinning of do- 
main walls, which pins the polarization in a particular 
direction, ra ther than any fundamental reduction of the 
polarisation. IScott and Poulignvl |)1988|) demonstrated in 
KNO3 via Raman spectroscopy that only a very small 
part of the sample was converted from the ferroelectric 
to non ferroelectric phase with fatigue, thus implying 
that fatigue must be caused by pinning of the domain 
walls. They also demonstrated that the domain walls 
could be de-pinned via the application of a large field, as 
shown in Fig. 17. The pinning of domain walls has also 
been observed directly by Ato mic Force microsco py by 
iGruverman et all l)l996|) and bv lColla et all l)l998ft . 

There is a fairly large body of evidence that oxygen va- 
cancies pla^_s^mekw_2art in the fatigue process. Auger 
data of lScott et all ()1991[) show areas of low oxygen con- 
centration in a region near to the metal electrodes, imply- 
ing a region of greater oxygen v acancy da t a. Scott et a l 
also reproduced Auger data from lTroeger I ^unpublished) 
for a film that had been fatigued by 10 10 cycles show- 
ing an increase in the width of the region with depleted 
oxygen near the platinum electrode (Fig. ITT)|) . There is 
however no corresponding change at the gold electrode. 
Gold does not form oxides, this might be an explana- 
tion of the different behaviour at the two electrodes. Al- 
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though some researchers believe that platinum also does 
not form oxides, the adsorption of oxygen onto Pt sur- 
faces is actually a large area of research, because of the 
important role platinum plays as a catalyst in fuel cell 
electrodes. Oxygen is not normally adsorbed onto gold 
surfac es but can be if there is significant surface rough- 
ness (|rvlills et all . \200$\ . 




FIG. 19 \Scott et alJ 1 199 A) fa) Auger depth profile of PZT 
thin film capacitor (b) Effect of fatigue on oxygen concentra- 
tion near the electrode 

It has also been found experimentally that films fa- 
tigue differently in atmosph eres conta i ning d if ferent oxy- 
gen p artial pressures. (jBrazier et all . fl999^ . iPan et al\ 
(1996) claim to have seen oxygen actually leaving a fer- 
roelectric sample d uring switching, though we note that 
iNuffer et all {20011 clai m this to be an ex perimental arte- 
fact. The results of ISchloss et all l)2002f) are very inter- 
esting in that they show directly by O 18 tracer studies 
that the oxygen vacancies redistribute themselves during 
voltage cycling. In their original paper they concluded 
that the redistribution of oxygen vacancies this is not 
the cause of fatigue, because they did not see redistribu- 
tion of O 18 when the sample has been annealed, though 
the sample still fatigues, however in a more recent pub- 
lication they conclude the reason they could not see the 
oxygen tracer distribution was more probably due to a 
change in the oxygen pe r meab ility of the electrode after 
annealing l|Schloss et qj].l2004h. 

It has been known for some time that the fatigue of 
PZT films can be improved by the use of oxi de electrodes, 
such a s iridium oxide or ruthenium oxides. lArauio et all 
(1995) explain the improved fatigue resistance by the fact 
that oxides of iridium and platinum can reduce or re- 
oxidise reversibly and repeatedly without degradation. 
It is for the same reason that iridium is preferred to plat- 
inum an a electrode materials for medical ap plications 
where this property was originally studied by iRobbleel 
(1986). This property does make the leakage current 
properties of these electrodes more complicated, and gen- 
erally films with Ir/Ir02 or Ru/Ru02 electrodes have 
higher leakage currents than those with platinum elec- 
trodes. Unless carefully annealed at a certain temper- 
ature RuC>2 electrodes will have elemental Ru metallic 
islands. Since the work function for Ru is 4.65 eV and 
that for Ru02 is 4.95 eV, almost all the current will pass 
through the Ru islands, producing h ot spots and occa- 
sional shorts l|Hartmann et all . |2000[) . By contrast, al- 
though one expects that there will be similar issues with 



mixtures of Ir and Ir02 in iridium based electrodes, when 
metallic Ir is oxidized to IrQ 2 its w orkfunction decreases 
to 4.23 eV (|Chalamala et aJll999l) . 

The idea that planes of oxygen vacancies perpendicu- 
lar to the polarisat ion direct i on co uld pin domain walls 
is originally due to iBrennanl |)1993|) . Subsequently, in a 
theoretical microscopic stud y of ox ygen vacancy defects 
in PbTiC>3 iPark and Chadil l|1998l) showed that planes 
of vacancies are much more effective at pinning do- 
main walls than single vac ancies. In bulk ferroelectrics 
lAr¥ and Neumann! I|1988|) have discussed how under 
repetitive cycling the vacancies can move from their 
originally randomly distributed sites in the perovskite 
structure to sites in planes parallel to the ferroelectric- 
electrode interface. We suspect that while this may ac- 
count for fatigue in bulk ferro electrics it is not th e oper- 
ative mechanism in thin films. IScott and Dawberl (2000) 
have suggested that in thin films the vacancies can reach 
sufficiently high concentrations that they order them- 
selves into planes in a similar way as occurs in Fe-doped 
bulk samples and on the surfaces of highly reduced spec- 
imens. Direct evidence that this occurs in bulk PZT was 
found usi ng Atomic force m i crosco pe imagery of PZT 
grains bv ILupascu and Rabd l)2002|) (Fig. Recently 
evidence of oxygen vacancy ordering has also been found 
in barium titanate redu ced after an accelerated life test 
l|Woodward et ql.U2004|) . 




FIG. 20 Atomic force microscopy images (Lupascu) (a) before 
and (b) after, cycling showing evidence of planes of oxygen 
vacancies in the fatigued sample. Image width = 10 /im. 

While most researchers acknowledge that oxygen va- 
cancies play a r o le in fatigue, it should be noted that 
iTaeantsev et all l)200l|) has aggressively championed a 
model of charge injection; however since this model is 
not developed into a quantitative form it is very hard to 
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verify or falsify it. Charge injection probably does play a 
role in fatigue, an idea at leas t in part supported b y the 
detailed experi mental studv oflDu and Cher] l(l998l) . but 
in the model of lDawb er and S cottl (l2000l> . f which draws 
upon the basic idea of hfocHmd""Desul 1)1992^) that fatigue 
is due to the electromigration of oxygen vacancies) it is 
not included, nevertheless most of the experimental re- 
s ults in the literature may be accounted for. The model 
of lDawber and Scottl l)200Cfl basically shows that in a fer- 
roelectric thin film under an AC field there is in fact a net 
migration of vacancies towards the interface and it is the 
high concentration of vacancies in this region that results 
in ordering of the vacancies and pinning of domain walls. 
The interfacial n ature of fat i gue in thin film s has been 
demonstrated bv lColla et all l|l998l) and bv lDimos et all 

To really understand fatigue better what is needed is 
more experiments that try to look at the problem in novel 
ways, standard electrical measurements alone probably 
cannot shed a great deal of additional light on the prob- 
lem, especially given the problems between comparing 
samples grown in different labs using different techniques. 
Very recently a very interesting study was undertaken by 
iDo et al\ l|2004j) using X-Ray microdiffraction to study 
fatigue in PZT. Using this technique they were able to 
see how regions of the film stopped switching as it fa- 
tigued. One of the key findings of this study was that 
there appears to be two fatigue effects operative, a re- 
versible effect that occurs when low to moderate fields 
are used for switching and an irreversible effect which 
occurs under very high fields. 



3. Retention failure 

Clearly a non-volatile memory that fails to retain the 
information stored in it will not be a successful device. 
Furthermore producers of memories need to be able to 
guarantee retention times over much longer periods of 
time then they can possibly test. A better understand- 
ing of retention failure is thus required so that models 
can used that allow a ccelerated r e tentio n tests to be car- 
ried out, the work of iKim et all fcODlf) is a step in this 
direction. 

It seems that imprint and retention failure are closely 
linked phenomena, ie. if a potential builds up across 
the system of the time it can destabilize the ferroelec- 
tric polarisation state and thus cause loss of information. 
Further comparison of retention-time data and fatigue 
data suggest quite strongly a link between the two ef- 
fects. DC degrada tion of resistance in BST seems also to 
be a related effect (Zaf ar et oil 119991. The electromigra- 
tion of oxygen vacancies under an applied field in a Fe- 
doped SrTiC-3 single cryst al has been directly observed 
via electro-chromic effects ijWaser et all . Il990|) . Oxygen 
vacancy redistribution under applied field has also been 
invoked to explain a slow relaxation of the capacitance 
in BST thin films ijBoikov et adl200lj) . It would seem to 



make sense that whereas fatigue relates to the cumulative 
motion of oxygen vacancies under an AC field, resistance 
degradation is a result of their migration under an ap- 
plied DC field and retention failure is a result of their 
migration under the depolarisation field or other built-in 
fields in the material. 



IV. FIRST PRINCIPLES 

With continuing advances in algorithms and computer 
hardware, first principles studies of the electronic struc- 
ture and structural energetics of complex oxides can now 
produce accurate, material-specific information relevant 
to the properties of thin-film ferroelectrics. In this sec- 
tion, we focus on first-principles studies that identify 
and analyze the characteristic effects specific to thin 
films. First, we briefly review the relevant methodologi- 
cal progress and the application of these methods to bulk 
ferroelectric materials. Next, we will survey the first- 
principles investigations of ferroelectric thin films and 
superlattices reported in the literature. It will be seen 
that the scale of systems that can be studied directly by 
first principles methods is severely limited at present by 
practical considerations. This can be circumvented by 
the construction of nonempirical models with parame- 
ters determined by fitting to the results of selected first- 
principles calculations. These models can be parame- 
terized interatomic potentials, permitting molecular dy- 
namics studies of nonzero temperature effects, or first- 
principles effective Hamiltonians for appropriate degrees 
of freedom (usually local polarization and strain). The 
form of the latter strongly resembles that of a Landau- 
Devonshire theory, providing a connection between first- 
principles approaches and the extensive literature on phe- 
nomenological models for the behavior of thin film ferro- 
electrics. The advantages and disadvantages of using first 
principles results rather than experimental data to con- 
struct models will be considered. In addition to allowing 
the study of systems far more complex than those that 
can be considered by first principles alone, this modelling 
approach yields physical insight into the essential differ- 
ences between bulk and thin film behavior, which will be 
discussed at greater length in subsection If V . B . 51 Finally, 
it will be seen that despite practical limitations, the com- 
plexity of the systems for which accurate calculations can 
be undertaken has steadily increased in recent years, to 
the point where films of several lattice constants in thick- 
ness can be considered. While this is still far thinner than 
the films of current technological interest, concommitent 
improvements in thin film synthesis and characterizaton 
have made it possible to achieve a high degree of atomic 
perfection in comparable ultrathin films in the context of 
research. This progress has led to a true relevance of cal- 
culational results to experimental observations, opening a 
meaningful experimental-theoretical dialogue. However, 
this progress in some ways only serves to highlight the 
full complexity of the physics of real ferroelectric films: as 
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questions get answered, more questions, especially about 
phenomena at longer length scales and about dynamics, 
are put forth. These challenges will be discussed in sub- 
section [EEO 



A. DFT studies of bulk ferroelectrics 

In parallel with advances in laboratory synthesis, the 
past decade has seen a revolution in the atomic-scale 
theoretical understanding of ferroelectricity, especially 
in perovskite oxides, through first-principles density- 
functional-theory (DFT) investigations. The central re- 
sult of a DFT calculation is the ground state energy 
computed within the Born-Oppenheimer approximation; 
from this the predicted ground state crystal structure, 
phonon dispersion relations, and elastic constants are 
directly accessible. The latter two quantities can be 
obtained by finite-difference calculations, or more effi- 
ciently, through the direct calculation of derivatives of 
the total energy through density-fun ctional perturbation 
theory (DFPT) l|Baroni et all 1200 ll) . 

For the physics of ferroelectrics, the electric polar- 
ization and its derivatives, such as the Born effective 
charges and the dielectric and piezoelectric tensors, are 
as central as the structural energetics, yet proper for- 
mulation in a first-principles context long proved to be 
quite elusive. Expressions for derivatives of the polar- 
ization corresponding to physically measurable quanti- 
ties were prese nted and applied in DFPT calculations in 
the late 1980's llde Gironcoli et aLlll989j) . A key concep- 
tual advance was establishing the correct definition of 
the electric polarization as a bulk property through the 
Berry phase formalism of King-Smith, Vanderbilt an d 
Resta (|King-Smith and Vanderbiltl Il993t iRestaL Il994|) - 
With this and the related Wannie r function expression 
ijKing-Smith and Vanderbiltl 1 1994|) . the spontaneous po- 
larization and its derivatives can be computed in a post- 
processing phase of a conventional total-energy calcula- 
tion, greatly facilitating studies of polarization-related 
properties. 

For perovskite oxides, the presence of oxygen and 
first row transition metals significantly increases the 
computational demands of density functional total en- 
ergy calculations compared to those for typical semi- 
conductors. Calculations for perovskite oxides have 
been reported using essentially all of the available first- 
principles methods for accurate representation of the 
electronic wavef unctions: all-electron methods, mainly 
linearized augmented plane wave (LAPW) and full- 
potential linearized augmented plane wave (FLAPW), 
linear muffin-tin orbitals (LMTO), norm-conserving 
and ultrasoft pscudopotentials, and projector-augmented 
wavefunction (PAW) potentials. The effects of dif- 
ferent choices for the approximate density functional 
have been examined; while most calculations are car- 
ried out with the local-density approximation (LDA), 
for many systems the effects of the generalized-gradient 



approximation ( GGA) an d weig hted-density approx- 
imation fWDAl llWu et all \ 2004) have been investi- 
gated, as well as the alternative use of the Hartree- 
Fock approach. Most calculations being currently 
reported are performed wit h an appropriate stan- 
dard package, mainly VASPllKresse and Furthmullerl 
1996; iKresse and Hafneil Il993j) ! with ultrasoft pseu - 
dopotentials and PAW potentials, ABINIT ijABJ l. 
with norm- conse rving pseudop otentials and PAW po- 
tentials, PWscf ijBaroni et ali) . with nor m-conserving 
and ultrasoft pseudopotentials, SIESTA ijSoler et all 
I2002D . with norm-conserving p seudopotentials, WIEN97 
CFLA PWl llBlaha et ^.1. 119901) and CRYSTAL (Hartree- 
Fock) l|Saunders et a/.l.ll999fl . 

To predict ground state crystal structures, the usual 
method is to minimize the total energy with respect to 
free structural parameters in a chosen space group, in 
a spirit similar to that of a Rietveld refinement in an 
experimental structural determination. The space group 
is usually implicitly specified by a starting guess for the 
structure. For efficient optimization, the calculation of 
forces on atoms and stresses on the unit cell is essential 
and is by now included in every standard first-principles 
implement ation following t he formalism of Hcllmann and 
Feynman (Feynman, 1939) fo r the forces and Nielsen and 
Martin ijNielsen and Martini Il985j) for stresses. 

The accuracy of DFT for predicting the 
ground state structures of ferroelectrics was 

first investigated for the prototypical ca ses of 

BaTi0 3 and P bTiO^ llCohen and Krakauerl Il990t 
ICohenl 1 19921 ICohen and Krakaueit Il992j) . and 
then exte nded to a larger class o f fer roelectric 
perovskites l|King-Smith and Vanderbiltl Il994|) . Ex- 
tensive studies of the structures of perovskite oxides 
and related ferro electric oxid e structures have since 
been carried out ijRestal l2003|) . The predictive power 
of first-principles calculation s is well illu strated by the 
results of Singh for PbZr0 3 ijSingbl F1995) . in which the 
correct energy ordering between ferroelectric and anti- 
ferroelectric structures was obtained and furthermore, 
comparisons of the total energy resolved an ambiguity 
in the reported space group and provided an accurate 
determination of the oxygen positions. 

It is important to note, however, that there seem to 
be limitations to the accuracy to which structural pa- 
rameters, particularly lattice constants, can be obtained. 
Most obvious is the underestimate of the lattice con- 
stants within the LDA, typically by about 1% (GGA 
tends to shift lattice constants upward, sometimes sub- 
stantially overcorrecting) . Considering that the calcula- 
tion involves no empirical input whatsoever, an error as 
small as 1% could be regarded not as a failure, but as a 
success of the method. Moreover, the fact that the under- 
estimate varies little from compound to compound means 
that the relative lattice constants and thus the type of lat- 
tice mismatch (tensile/compressive) betwen two materi- 
als in a heterostructure is generally correctly reproduced 
when using computed lattice parameters. However, for 
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certain questions, even a 1% underestimate can be prob- 
lematic. The ferroelectric instability in the perovskite 
oxides, i n particular, is known to be very sensitive to 
pressure l|Samaralll987|) and thus to the lattice constant, 
so that 1% can have a significant effect on the ferroelec- 
tric instability. In addition, full optimization of all struc- 
tural parameters in a low symmetry space group can in 
some cases, PbTi03 being the most well studied example, 
lead to an apparently spurious prediction, though fixing 
the lattice constants to their experimental values leads 
to good agreement for th e other structural parameters 
iSagh i-Szabo et q^J . ll998j) . Thus, it has become accept- 
able, at least in certain first-principles contexts, to fix the 
lattice parameters or at least the volume of the unit cell, 
to the experimental value, when this value is known. 

In a first-principles structural prediction, the initial 
choice of space group may appear to limit the chance 
that a true ground state structure will be found. In gen- 
eral, once a minimum is found, it can be proved (or not) 
to be a local minimum by computation of the full phonon 
dispersion and of the coupling, if allowed by symmetry, 
between zone-center phonons and homogeneous strain 
llGarcia and VanderbiHi Il99^. Of course, this does not 
rule out the possibility of a different local minimum with 
lower energy with an unrelated structure. 

For ferroelectrics, the soft-mode theory of ferroelec- 
tricity provides a natural conceptual framework for the 
identification of low-symmetry ground state structures 
and for estimating response functions. The starting point 
is the identification of a high-symmetry reference struc- 
ture. For perovskite compounds, this is clearly the cubic 
perovskite structure, and for layered perovskites, it is 
the nonpolar tetragonal structure. The lattice instabil- 
ities of the reference structure can be readily identified 
by the fir st-principles calculation of phonon dispersion 
relations llGhosez et all Il999t ISai and Vanderbiltt l2000t 
IWaehmare and Rabel Il997aft this being especially effi- 
cient within density functional perturbation theory. In 
simple ferroelectric perovskites, the ground state struc- 
ture is obtained to lowest order by freezing in the most 
unstable mode (a zone-center polar mode). This pic- 
ture can still be useful for more complex ground state 
structures that involve the freezin g in of two or more 
coup le d modes (e.g. PbZ r CK) l|Cockavne and Rabel 
l2000i IWaghmare and Rabel Il997bft . as well as for 
identifying low-energy structures th at might be sta- 
bilized by temperatu re or pressure l)Fennie and Rabel 
IStachiotti et all 12000ft The polarization induced by the 
soft polar mode can be obtained by computation of 
the Born effective charges, yielding the mode effective 
charge. The temperature-dependent frequency of the 
soft polar mode and its coupling to strain are expected 
largely to determine the dielectric and piezoelectric re- 
sponse of ferroelectric and near-ferroelectric perovskite 
oxide s; this idea has been the basis of several calcula- 
tions llCockavne and Rabel H^ lc arcia and Vanderbihl 
11998ft . 

Minimization of the total energy can similarly be 
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non-stoichiometric planar defects such 
as antiphase doma i n bo undaries l)Li et all 120021 
ISuzuki and Fuiimotol l200lft .in bulk perovskite oxides. 
The supercell used must accommodate the defect ge- 
ometry, and generally must contain many bulk primi- 
tive cells to minimize the interaction of a defect with its 
periodically-repeated images. Thus, these calculations 
are extremely computationally intensive, and many im- 
portant questions remain to be addressed. 

While much of the essential physics of ferroelectrics 
arises from the structural energetics, the polarization and 
the coupling between them, there has been increasing 
interest in ferroelectric oxides as electronic and optical 
materials, for which accurate calculations of the gap and 
dipole matrix elements are important. Furthermore, as 
we will discuss in detail below, the bandstructures enter 
in an essential way in understanding the charge trans- 
fer and dipole layer formation of heterostructures involv- 
ing ferroelectrics, other insulators, and metals. While 
density functional theory provides a rigorous foundation 
only for the computation of Born-Oppenhcimer ground- 
state total energies and electronic charge densities, it 
is also often used for investigation of electronic struc- 
ture. In the vast majority of density functional im- 
plementations, calculation of the ground state total en- 
ergy and charge density involves the computation of a 
bandstructure for independent electron states in an ef- 
fective potential, f ollowi ng the work of Kohn and Sham 
ijKohn and Shaml [l965). This bandstructure is gener- 
ally regarded as a useful guide to the electronic structure 
of mate rials, including perovskite and l ayered perovskite 
oxides ijCohenl Il992t iRobertson et all Il996t iTsai et all 
12003ft . It should be noted that it is consistently found 
that with approximate functionals, such as the LDA, 
the fundamental bandgaps of insulators and semiconduc- 
tors, including perovskite ferroelectrics, are substantially 
underestimated. While for narrow gap materials, the 
system may even be erroneously found to be metallic, 
for wider gap systems such as most of the simple ferro- 
electric perovskite compounds considered here, the band 
gap is still nonzero and thus the structural energetics in 
the vicinity of the ground state structure is unaffected. 
While this error might be considered to be an insuperable 
stumbling block to first-principles investigation of elec- 
tronic structure and related properties, there is at present 
no widely available, computationally tractable, alterna- 
tive (there are, though, some indications that the use of 
exa ct-exchange functionals can eliminate much of this er- 
ror ()Piskunov et all\2004) : it has long been known that 
Hartree-Fock, i.e. exact exchange only, leads to overesti- 
mates of the gaps). The truth is that, as will be discussed 
further below, these results can, with care, awareness of 
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the possible limitations and judicious use of experimen- 
tal input, be used to extract useful information about the 
electronic structure and related properties in individual 
material systems. 

At present, the computational limitations of full first- 
principles calculations to 70-100 atoms per supercell have 
stimulated considerable interest in the development and 
use in simulations of effective models, from which a 
subset of the degrees of freedom have been integrated 
out. Interatomic shell-model potentials have been de- 
veloped for a number of perovskite oxide systems, elim- 
inating most of the electronic degrees o f freedom ex- 
cept for those represent e d by the shells llHeifets et all 
l2000bt ISepliarskv eTdl 12004 iTinte et all Il999jh A 
more dramatic reduction in the number of degrees of 
freedom is performed to obtain effective Hamiltonians, 
in which typically one vector degree of freedom decribes 
the local polar distortion in each unit cell. This ap- 
proach has proved useful for describing finite temperature 
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extent that the parameters appearing in these potentials 
are determined by fitting to selected first-principles re- 
sults (structures, elastic constants, phonons), these ap- 
proaches can be regarded as approximate first-principles 
methods. They allow computation of the polarization as 
well as of the structural energetics, but not, however, of 
the electronic states. Most of the effort has been focused 
on BaTi0 3 , though other perovskites, including PbTi0 3 
and KNb0 3 , SrTi0 3 and (Ba,Sr)Ti0 3 and Pb(Zr,Ti)0 3 
have been investigated in this way, and useful results ob- 
tained. 



B. First-principles investigation of ferroelectric thin films 



havior can best be achieved by detailed quantitative ex- 
amination of individual systems combined with the con- 
struction of models incorporating various aspects of the 
physics, from which more general organizing principles 
can be identified. In first principles calculations, there 
is a freedom to impose constraints on structural param- 
eters and consider hypothetical structures that goes far 
beyond anything possible in a real system being studied 
experimentally. This will allow us to isolate and examine 
various influences on the state of a thin film: epitaxial 
strain, macroscopic electric fields, surfaces and interfaces, 
characteristic defects associated with thin-film growth, 
and "true" finite size effects, and how they change the 
changes in atomic arrangements, electronic structure, po- 
larization, vibrational properties, and responses to ap- 
plied fields and stresses. While the main focus of this 
review is on thin films, this approach also applies nat- 
urally to multilayers and superlattices. Extending our 
discussion to include these latter systems will allow us to 
consider the effects of the influencing factors in different 
combinations, for example the changing density of inter- 
faces, the degree of mismatch strain, and the polarization 
mismatch. These ideas also are relevant to investigating 
the behavior of bulk layered ferroelectrics, which can be 
regarded as natural short-period superlattices. 

Within this first-principles modelling framework, we 
can more clearly identify specific issues and results for 
investigation and analysis. The focus on modelling is 
also key to the connection of first-principles results to the 
extensive literature on phcnomenological analysis and to 
experimental observations. This makes the most effective 
use of first-principles calculations in developing a concep- 
tual and quantitative understanding of characteristic thin 
film properties, as manifested by thickness dependence as 
well as by the dependence on choice of materials for the 
film, substrate and electrodes. 



The fascination of ferroelectric thin films and superlat- 
tices lies in the fact that the properties of the system as a 
whole can be so different from the individual properties 
of the constituent material(s). Empirically, it have been 
observed that certain desirable bulk properties, such as 
a high dielectric response, can be degraded in thin films, 
while in other investigations there are signs of novel in- 
teresting behavior obtained only in thin film form. 

Theoretical analysis of the observed properties of thin 
films presents a daunting challenge. It is well known that 
the process of thin film growth itself can lead to non- 
trivial differences from the bulk material, as observed 
in studies of homoep itaxial oxide films such as SrTi0 3 
ijKlenov et al .], 12004 . However, as synthetic methods 
have developed, the goal of growing nearly ideal, atom- 
ically ordered, single crystal films and superlattices is 
coming within reach, and the relevance of first-principles 
results for perfect single crystal films to experimental ob- 
servations emerging. 

As will be clear from the discussion in the rest of this 
section, an understanding of characteristic thin film be- 



1. First principles methodology for thin films 

The fundamental geometry for the study of thin films, 
surfaces and interfaces is that of an infinite single- 
crystalline planar slab. Since three-dimensional periodic- 
ity is required by most first-principles implementations, 
in those cases the slab is periodically repeated to pro- 
duce a supercell; a few studies have been carried out 
with electronic wavefunction basis sets that permit the 
study of an isolated slab. The variables to be specified 
include the orientation, the number of atomic layers in- 
cluded, choice of termination of the surfaces, and width 
of vacuum layer separating adjacent slabs. As in the 
first-principles prediction of bulk crystal structure, choice 
of a space group for the supercell is usually established 
by the initial structure; relaxations following forces and 
stresses do not break space group symmetries. The di- 
rection of the spontaneous polarization is constrained by 
the choice of space group, allowing comparison of unpo- 
larized (paraelectric) films with films polarized along the 
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normal or in the plane of the film. 

As we will see below, in most cases the slabs are very 
thin (ten atomic layers or fewer). It is possible to re- 
late the results to the surface of a semi-infinite system 
or a coherent epitaxial film on a semi-infinite substrate 
by imposing certain constraints on the structures con- 
sidered. In the former case, the atomic positions for 
interior layers are fixed to correspond to the bulk crys- 
tal structure. In the latter, the in-plane lattice parame- 
ters of the supercell are fixed to the corresponding bulk 
lattice parameters of the substrate, which is not other- 
wise explicitly included in the calculation. More sophis- 
ticated methods developed to deal with the coupling of 
vibration al modes at the su r face w ith bulk modes of the 
substrate ljLewis and Rappel Il996|) could also be applied 
to ferroelectric thin films, though this has not yet been 
done. 

The LDA underestimate of the equilibrium atomic vol- 
ume will in general also affect slab calculations, and sim- 
ilar concerns arise about the coupling of strain and the 
ferroelectric instability. As in bulk crystal structure pre- 
diction, it may in some cases be appropriate to fix certain 
structural parameters according to experimental or bulk 
information. In the case of superlattices and supercells 
of films on substrates, it may on the other hand be a 
good choice to work consistently at the (compressed) the- 
oretical lattice constant, since the generic underestimate 
of the atomic volume ensures that the lattice mismatch 
and relative tensile/compressive strain will be correctly 
reproduced. This applies for example to the technique 
mentioned in the previous paragraph, in which the ef- 
fects of epitaxial strain are investigated by performing 
slab calculations with an appropriate constraint on the 
in-plane lattice parameters. 

As in first-principles predictions of bulk crystal struc- 
ture, the initial choice of space group constrains, to a 
large extent, the final "ground-state" structure. If the 
supercell is constructed by choosing a bulk termination, 
the energy minimization based on forces and stresses will 
preserve the initial symmetry, yielding information about 
surface relaxations of the unreconstructed surface. A 
lower-energy structure might result from breaking addi- 
tional point or translational symmetries to obtain a sur- 
face reconstruction. This type of surface reconstruction 
could be detected by computing the Hessian matrix (cou- 
pled phonon dispersion and homogeneous strain) for the 
relaxed surface. More complex reconstructions involv- 
ing adatoms, vacancies or both would have to be stud- 
ied using appropriate starting supercells. Information re- 
garding the existence and nature of such reconstructions 
might be drawn from experiments and/or from known 
reconstructions in related materials. 

One very important consideration in the theoretical 
prediction of stable surface orientations and termina- 
tions and of favorable surface reconstructions is that 
these depend on the relative chemical potential of the 
constituents. Fortunately, since the chemical poten- 
tial couples only to the stoichiometry, the prediction of 



the change of relative stability with chemical potential 
can be made with a single total-e nergy calculation f or 
each structure (see, for example, l|Mever et all H999)). 
Because of the variation in stoichiometry for different 
(001) surface terminations, what is generally reported is 
the average surface energy of symmetric AO and BO2- 
tcrminated slabs. 

A problem peculiar to the study of periodically re- 
peated slabs with polarization along the normal is the 
appearance of electric fields in the vacuum. As shown in 
Fig. 1211 this occurs because there is a nonzero macro- 
scopic depolarizing field in the slab and thus a nonzero 
potential drop between the two surfaces of the slab. As 
the potential drop across the entire supercell must be 
zero, this inevitably leads to a nonzero electric field in the 
vacuum. Physically, this can be interpreted as an exter- 
nal field applied uniformly across the unit cell which acts 
partially to compensate the depolarizing field. An anal- 
ogous situation arises for an asymmetrically terminated 
slab when the two surfaces have different work functions. 
To eliminate the artificial field in the vacuum, one tech- 
nique is to introduce a dipole l ayer in the mid-v acuum 
region far away from the slab l)Bengtssonlll999|) . This 
can accommodate a potential drop up to a critical value, 
at which point electrons begin to accummulate in an ar- 
tificial well in the vacuum region (see Fig. 12 ip . This 
approach can be also be used to compensate the depo- 
larization field in a perpendicular polarized film, though 
it may happen that the maximum field that can be ap- 
plied is smaller than that needed for full compensation. 
Alternatively, by using a first-principles implementation 
with a localized basis set, it is possible to perform com- 
putations for isolated slabs and thus avoid not only the 
spurious electric fields, but also the interaction between 
periodic slab images present even for symmetric nonpo- 
lar slabs. Comparison between results obtained with the 
two approaches is presented in l)Fu et all\l999\i . 

In determining the properties of an ultrathin film, 
the film-substrate and/or film electrode interface plays 
a role at least as important as the free surface. In first- 
principles calculations, the atomic and electronic struc- 
ture of the relevant interface(s) is most readily obtained 
by considering a periodic multilayer geometry identical 
to that used for computing the structure and properties 
of superlattices. To simulate a semi-infinite substrate, 
the in-plane lattice constant should be fixed to that of 
the substrate bulk. The relaxation of a large number of 
structural degrees of freedom requires substantial com- 
puter resources, and some strategies for efficiently gener- 
ating a good starting structure will be described in the 
discussion of specific systems in Scction liV.B.3l 

Calculations of quantities characterizing the electronic 
structure are based on use of the Kohn-Sham one- 
electron energies and wave functions. Band structures for 
lxl (001) slabs are generally displayed in the 2D surface 
Brillouin zone for k superce u -z =0. One way to identify 
surface states is by comparison with the appropriately 
folded-in bulk bandstructure. Another analysis method 
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FIG. 21 Schematic picture of the planar-averaged potential 
v(z) for periodically repeated slabs: (a) with periodic bound- 
ary conditions, (b) potential of the dipole layer, (c) dipole- 
corrected slabs with vanishing external electric field, and (d) 
dipol e-corrected slabs with vanish ing internal electric field. 
From lMever and Vandcr biltL l200ll 



is to compute the partial density of states projected onto 
each atom in the slab. From these plots, the domi- 
nant character of a state at a particular energy can be 
found. In addition, an estimate of valence band offsets 
and Schottky barriers at an interface can be obtained 
by analyzing the partial density of states for a superlat- 
tice of the two constituents. This is done comparing the 
energies of the highest occupied states in the interior of 
the relevant constituent layers (because of the band gap 
problem, the positions of the conduction bands are com- 
puted by using the experimental bulk bandgaps) . This 
estim ate can be refined, as described in lJunauera et ali . 
l2003t by computation of the average electrostatic energy 
difference between the relevant constituent layers. 

As the computational resources required for full first- 
principles calculations even of the simplest slab-vacuum 
system are considerable, there is a strong motivation to 
turn to interatomic potentials and effective Hamiltoni- 
ans. Interatomic potentials based on shell models fit- 
ted to bulk structural energetics are generally directly 
transferred to the isolated slab geometry, with no changes 
for the undercoordinated atoms at the surface. As will 
be discussed further below, this approach seems to be 
successful in reproducing the relaxations observed in full 
first-principles studies and has been applied to far larger 
supercells and superlattices. In the case of effective 



Hamiltonians, it is, at least formally, possible simply to 
perform a simulation by removing unit cells and use bulk 
interaction parameters for the unit cells in the film. For 
a more accurate description, modification of the effective 
Hamiltonian parameters for the surface layers is advis- 
able to restore the charge neutrality sum rule for the 
film llGhosez and Rabel 120001 iRuini et a/.LIl998fl . In ad- 
dition, the effects of surface relaxation would also result 
in modified interactions at the surface. 



2. Overview of systems 

In this section, we present a list of the materials and 
configurations that have been studied, followed by a brief 
overview of the quantities and properties that have been 
calculated in one or more of the reported studies. A more 
detailed description of the work on individual systems is 
provided in the following subsection. 

The configurations that have been considered to 
date in first principles studies can be organized into 
several classes. The simplest configuration is a slab 
of ferroelectric material alternating with vacuum; 
this can be used to investigate the free surface of 
a semi-infinite crystal, an unconstrained thin film, 
or an epitaxial thin film constrained to match the 
lattice constant of an implicit substrate. Specific 
materials considered included BaTiOs ( (001) surface s 
ilCohenl Il99fl Il997t IFu et all Il999t iHeifets et al 




2001 
2001 



L 
it. 



iKr^rnar and Fd [2003: iMevertmc^aiiTterbii' 
an d (110) su rfa ces IIHeifets et all 1 2001a) 
SrTiO^ llHeifets et all l2002al l2001al 12002 
iKubo and NozoveL 12003 iPadilla and Vanderbihl I199S.. 
PbTiQ 3 Emigarogm^r lab^r iGn^s^^m^^aba I200C 
Mever and Vanderbiltl 120011) KNb0 3 (|Heifets et all 
20004 and KTaCK l)Li et all l2003a|l Another type 
of configuration of comparable complexity is ob- 
tained by replacing the vacuum by a second material. 
If this is another insulating perovskite oxide, the 
calculation can yield information a bout ferroelectric 
dielectric (e.g BaTiCh/SrTi CK llNeaton and Rabe, 



2003) and KNb0 3 /KTa0 3 llSepliarskv et all 



2002, 



2001)) or ferroelectric- ferroelectric interfaces and su- 
perlattices. This configuration can also be used to 
study the interface between the ferroelectric and a 
dielectric (nonferro electric) oxide ( e.g. BaTi0 3 /BaO 
and SrTi0 3 /SrO ()Junauera et all l2003|) ). Replace- 
ment of the vacuum by a conductor simulates a film 
with symmetric al top and bottom e lectro des, e.g. 
BaTiO^/SrRuO ^Junauera and Ghosei l2003j) . More 
complex multilayer geometries including two or more 
different materials as well as vacuum layers have been 
used to simulate ferroelectric thin film interactions with 
the substrate (e.g. PbTi0 3 /SrTi0 3 /PbTi0 3 /vacuum 
ll Johnston and Rabel) . and wi th realisti c elec trodes 
(e.g. Pt/BaTiO^/Pt/ vacuumllRao et all Il997|) and 
Pt/PbTi0 3 /Pt /vacuum (|Sai and RarjpeD . as well as the 
structure of epitaxial alkaline-earth oxide on silicon, 
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Ba/Pb • Ti O O 

FIG. 22 Schematic illustr ation of the structure of the first 
three surface layers. From lMever and Vanderbihl l200ll . 

used as a buffer layer fo r growth of perovskite oxide 
films l)McKee et al\ . \200$i . 

In each class of configurations, there are corresponding 
quantities and properties that are generally calculated. 
For the single slab-vacuum configuration, for each ori- 
entation and surface termination the surface energy is 
obtained. While in some of the initial studies the atomic 
positions were fixed according to structural information 
from the bulk llCohen! [1996), in most current studies 
relaxations are obtained by energy minimization proce- 
dures. For the most commonly studied perovskite (001) 
slab, the relaxation geometry is characterized by changes 
in interplane spacings and rumplings quantified as shown 
m Fig. El Most studies assume bulk periodicity in the 
plane. For the study of surface reconstructions, it is nec- 
essary to expand the lateral unit cell, leading to a sub- 
stantial additional cost in computational resources. Most 
attention has been focused on the paraelectric SrTiC>3, 
although recently studies have been carried out as well 
for PbTi03. A byproduct of the total energy calcula- 
tion that is often though not universally presented is the 
bandstructure and or density of states; the local density 
of states at the surface is of particular interest. 

The two-component superlattice configuration can be 
taken to model a film on a substrate and/or with elec- 
trode layers, or to model an actual superlattice such as 
that which can be obtained by techniques such as MBE. 
In these studies, the interface and coupling between the 
two constituents is of primary interest; combinations 
most considered to date are ferroelectric + paraelectric, 
ferroelectric + dielectric, and ferroelectric + metal, while 
the combination of two ferroelectrics or ferroelectric + 
ferromagnetic material has been less intensively investi- 
gated. The main questions of interest are the structural 
rearrangements at the interface and the change in the 
structure, polarization, and related properties of individ- 
ual layers relative to the bulk resulting from the interac- 
tion with the other constituents. Analysis of the trends 
with varying thickness (es) of the ferroelectric film and, in 
the superlattice, other constituents is particularly useful. 
The electronic structure of these systems can be most 
readily characterized by the band offset between the two 
materials, which should also control the charge transfer 
across the interface, formation of a dipole layer, and the 



potential difference between the two constituents. In the 
case of a metal, this will determine the type of contact. 
The existence of interface states is also very relevant to 
the physical behavior of the system. 

At a considerable increase in computational expense 
but also in realism, a system with three or more compo- 
nents can be studied; e.g. the combination of a substrate, 
a film, and vacuum. The main questions of interest in the 
few such studies to date are the analysis of the ferroelec- 
tric instability in the film, and the film-induced changes 
in the substrate layers closest to the interface. As in two- 
component heterostructures, the partial density of states 
and the layer-average electrostatic potential are also use- 
ful in extracting the electronic behavior of the system. 

In all of these studies, one of the main questions is 
that of the structure and polarization of the ferroelec- 
tric layer compared to that of the bulk. Certainly, the 
change of environment (electrical and mechanical bound- 
ary conditions) and the finite dimensions (film thickness, 
particle size) are expect strongly to affect the structure 
and perhaps to eliminate the ferroelectric instability en- 
tirely. Relevant quantities to examine include the rela- 
tive stability of lower and higher symmetry phases, spa- 
tial variation in polarization, changes in the average po- 
larization magnitude and direction, and the depth and 
shape of the ferroelectric double-well potential. These 
changes can also be expected to lead to changes in the 
dielectric and piezoelectric response of thin films and su- 
perlattices, which can be studied theoretically and com- 
pared with experiments. The implications of the various 
first-principles studies included in this review will be de- 
scribed at further length below. 



3. Studies of individual one component-systems 

In this section, we describe a representative sampling of 
first principles studies and their results. Most of the lit- 
erature has concentrated on BaTiOs, providing a useful 
comparative test of various first-principles implementa- 
tions, as well as a benchmark for evaluation and analysis 
of results on other systems. We first consider the calcu- 
lations for single slabs of pure material that focus on the 
properties of surfaces: surface relaxation, surface recon- 
structions and surface states. Depending on the struc- 
tural constraints, these calculations are relevant cither 
to the surface of a semi-infinite bulk, for a free-standing 
thin film, or for a thin film epitaxially constrained by a 
substrate. This will be followed by discussions of studies 
of systems with two or more material components. 



a. BaTiOs For BaTiC>3, full first-principles results have 
been reported primarily for the (001) orientation, with a 
few results for the (110) and (111) orientations. In the 
slab- vacuum configuration, systems up to 10 atomic lay- 
ers have been considered. The unpolarized slab is com- 
pared with slabs with nonzero polarization, in the plane 
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and/or along the normal. After reviewing the results on 
structures, we will describe the results of extension to 
larger-scale systems through the use of interatomic po- 
tentials. The discussion of BaTiC>3 will be concluded by 
description of the first-principles results for surface elec- 
tronic structure. 

First-principles FLAPW calculations for the BaTiO 
1996), a nd lat er extended using the LAPW+LO method 



(001) and (111) slabs were first presented in 1995 



BaTiOs 
llCoherl 



iCoheillT997h . The supercells contained 6 and 7 atomic 
layers, corresponding to asymmetric termination and two 
symmetric terminations (BaO and Ti02), and an equal 
vacuum thickness. The central mirror plane symmetry 
z— > -z symmetry is broken for the asymmetric termina- 
tion even in the absence of ferroelectric distortion. The 
primitive cell lattice constant was fixed at the experimen- 
tal cubic value 4.01 A. Total energies of several selected 
paraclcctric and ferroelectric structures were computed 
and compared: the ideal paraelectric slab and ferroelec- 
tric slabs with displacements along z corresponding to 
the experimental tetragonal structure. For the asymmet- 
rically terminated slab, both the ferroelectric structure 
with polarization towards the BaO surface (+) and the 
other towards the Ti0 2 surface (-) were considered. The 
surface layers were relaxed for the ideal and (+) ferroelec- 
tric asymmetrically-terminated slabs and the ferroelec- 
tric BaO-BaO slab. It was found that the depolarization 
field of the ferroelectric slabs indeed strongly destabilizes 
the ferroelectric state, as expected, even taking into ac- 
count the energy lowering due to the surface relaxation. 
In all slabs considered, this consists of an inward-pointing 
dipole arising from the relative motion of surface cations 
and oxygens. The average surface energy of the ideal 
BaO and Ti0 2 surfaces is 0.0574 eV/A 2 =0.923 eV per 
surface unit cell. 

In iPadilla and Vanderbilti . Il997l ultrasoft pseudopo- 
tential calculations with fully relaxed atomic coordinates 
were reported for symmetrically terminated (both BaO 
and Ti02) 7-layer BaTiOs (001) slabs separated by two 
lattice constants of vacuum. The in-plane lattice con- 
stant was set equal to the theoretical equilibrium lat- 
tice constant a computed for the bulk tetragonal phase 
(a = 3.94 A). The average surface energy of the ideal 
BaO and Ti02 surfaces is 1.358 eV per sur face un i t cell; 
at least part of the difference relative to ICohenl Il996l 
could be due to the different lattice constant. Relax- 
ations were reported for unpolarized slabs and for polar- 
ized slabs with polarization along (100) (in the plane of 
the slab). Deviations from the bulk structure were con- 
fined to the first few atomic layers. The surface layer 
relaxes substantially inwards, and rumples such that the 
cation (Ba or Ti) moves inward relative to the oxygens, 
as in Fig. [221 While the relaxation energy was found to 
be much greater than the ferroelectric double well depth, 
the in-plane component of the unit-cell dipole moment 
was seen to be relatively insensitive to the surface relax- 
ation, with a modest enhancement at the Ti0 2 termi- 
nated surface and a small reduction at the BaO termi- 



nated surface. The relative stability of BaO and Ti02 
terminations were compared and both found to be stable 
depending on whether the growth was under Ba-rich or 
Ti-rich conditions. 

This investigati on was extended in 

iMever and Vanderbiltl l200ll to 7-layer and 9-layer 
polarized slabs with polarization along the normal. 
The problem of the artificial vacuum field in this 
periodically repeated slab calculation was addressed 
by the techniques of introducing an external dipole 
layer in the vacuum region of the supercell described 
above in Section IIV.B.1I This technique can also be 
used to generate an applied field that partially or fully 
compensates the depolarization field for BaTiOs slabs. 
As a function of applied field, the change in structure 
can be understood as arising from oppositely-directed 
electrostatic forces on the positively charged cations 
and negatively charged anions, leading to corresponding 
changes in the rumplings of the atomic layers and 
field-induced increases of the layer dipoles. Analysis of 
the internal electric field as a function of the applied 
field allows a determination of whether the slab is 
paraelectric or ferroelectric. The BaO-terminated slab 
is clearly ferroelectric, with vanishing internal electric 
field at an external field of 0.05 a.u. and a polarization 
of 22.9 \iCcmr" 1 , comparable to the bulk spontaneous 
polarization. The ferroelectric instability is suppressed 
in the Ti02-terminated slab, which appears to be 
marginally paraelectric. 

meth od wa s use d in 
and iFu et al\ l|l999|) . In 
a detailed analysis of the 
bonding was performed using tightbinding parametriza- 
tion. For the 7-layer BaO-terminated slab, the reported 
displacement of selected Ti an d O atoms is in good 
agreem ent with the results of IPadilla and Vanderbiltl 
and these calculations were extended to slabs of 
up to 15 layers. Fu et al. l)Fu et all Il999|) performed 
Hartree-Fock calculations for slabs of two to eight atomic 
layers, with symmetric and asymmetric terminations. 
Using a localized basis set, they were able to perform 
calculations for isolated slabs as well as periodically 
repeated slabs. Calculations of the macroscopically 
averaged planar charge density, surface energy and 
surface dynamical charges were reported as a function 
of thickness and termination for a cubic lattice constant 
of 4.006 A. The relative atomic positions were fixed 
to their bulk tetragonal structure values (note that 
this polarized structure in both isolated and periodic 
boundary conditions has a very high electrostatic energy 
and is not the ground state structure). This would 
significantly affect the comparison of the computed 
surface properties with experiment. In particular, it is 
presumably responsible for the high value of the average 
surface energy reported (1.69 eV per surface unit cell). 
However, a useful comparison between isolated and 
periodic slabs is possible. It was found that the surface 
charge and surface dipoles of isolated slabs converge 
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quite rapidly as a function of slab thickness and can be 
used, combined with a value of taken from the bulk, 
to extract a spontaneous polarization of 0.245 C/m 2 
(corrected to zero field using the electronic dielectric 
constant — 2.76), to be compared with 0.240 C/m 2 
from a Berry-phase calculation. This is only slightly 
less than the bulk value of 0.263 C/m 2 taken from 
experiment. The average surface energy for the two 
terminations of symmetrically terminated slabs is 0.85 
eV per surface unit cell. Surface longitudinal dynamical 
charges differ considerably from bulk values, satisfying 
a sum rule that the dynamical charges at the surface 
planes add up to half of the corresponding bulk value 
l|Ruini et a/l ll998). Convergence of all quantities with 
slab and vacuum thickness of periodically repeated slabs 
was found, in comparision, to be slow, with significant 
corrections due to the fictitious field in the vacuum (for 
polarized slabs) and the interaction between slab images. 

The isolated slab was also t he subject of an FLAPW 
study l|Krcmar and Fii l2003h . The symmetric Ti0 2 - 
Ti02 (9-layer) and asymmetric Ti02-BaO (10-layer) 
slabs were considered in a paraelectric structure with a 
fixed to 4.00 A and a polar tetragonal structure with 
a and c = 4.00 and 4.04 A, respectively. For the cu- 
bic Ti02-terminated slab, displacements in units of c 
for the surface Ti, surface O, subsurface Ba and sub- 
surface O are -0.021, +0.0 07, +0.022 and -0.009 c to be 
compared with the results ljPadilla and Vanderbiltt Il997|) 
(-0.0389,-0.0163,+0.0131,-0.0062) for a periodically re- 
peated 7-layer slab with lattice constant 3.94 A. The 
tetragonal phase was relaxed to a convergence criterion 
of 0.06 eV/A on the atomic forces; the rumplings of 
the layer s follow overall the same pa ttern as that re- 
ported in iMever and Vanderbiltl l200li with an inward- 
pointing surface dipole arising from surface relaxation, 
though the reduction of the rumpling in the interior is 
not as pronounced as for th e zero-applied field case in 
IMever and Vanderbiltl l200ll The energy difference be- 
tween the paraelectric and ferroelectric slabs was not re- 
ported. 

With interatomic potentials, it is possible to study 
additional aspects of surface behavior in BaTiOs thin 
films and nanocrystals. The most important feature of 
interatomic studies of thin films relative to full first- 
principles calculations is the relative ease of extending 
the supercell in the lateral direction, allowing the forma- 
tion of 180 degree domains and mol ecular dynamics stud- 
ies of finite temperature effects. In lTinte and Stachiottil 
l)200lj) . a 15-layer Ti02-terminated slab periodically re- 
peated with a vacuum region of 20 A was studied, using 
interatomic potentials that had previously been bench- 
marked against first-principles sur face relaxations and 
energies l|Tinte and Stachiottil l2000|) . The unconstrained 
(stress-free) slab is found to undergo a series of phase 
transitions with decreasing temperature, from a paraelec- 
tric phase to ferroelectric phases, first with polarization 
in the plane along (100), and then along (110). Enhance- 
ment of the surface polarization at low temperatures ap- 
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FIG. 23 Cell-by-cell out-of-plane (top panel) and in-plane 
(bottom panel) polarization profiles of a randomly chosen 
chain perpendicular to the slab surface for different misfit 
strains r\ at T = K. In the in-plane p olarization profiles 
p v —p x . From lTinte and Stachiottil 1200 ll . 



pears to be linked to the existence of an intermediate 
temperature regime of surface ferroelectricity. For slabs 
with a strongly compressive epitaxial strain constraint (?y 
= -2.5%), there is a transition to a ferroelectric state with 
180 degree domains and polarization along the surface 
normal. At the surface, the polarization has a nonzero x 
component and a reduced z-component, giving a rotation 
at the surface layer. The width of the stripe domains can- 
not be determined, as it is limited by the lateral supercell 
size at least up to 10x10. Reducing the compression to r\ 
= -1.0% also gives a 180 degree domain structure in the z 
component of the polarization, combined with a nonzero 
component along [110] in the interior of the film as well 
as the surface, as can be see from examination of Fig. 
1231 The thickness dependence of the transition temper- 
ature to this ferroelectric domain phase was studied at 
7] = -1.5%, with T c decreasing from the 15-layer film to 
the 13-layer and lowest T c 11-layer film. At and above a 
critical thickness of 3.6 nm, stress-free films exhibit the 

same ferroelectric-domain ground-state s t ructur e. 

Heifets and coworkers l)Heifets et all l2001al [2000b) 
studied BTO (001) and (110) surfaces of an i solated slab 
using the shell model of lHeifets et all l2000bl Between 1 
and 16 atomic planes were relaxed in the electrostatic 
potential of a a rigid slab of 20 atomic planes whose 
atoms were fixed in their perfect (presumably cubic) lat- 
tice sites. The relaxed structures of the two (001) termi- 
nations are in good agreement with other calculations, 
except for the sign of the surface dipole in the relax- 
ation of the BaO terminated surface, which is found to 
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be positive (though small). The (110) surfaces are found 
to have much higher surface energy, except for the re- 
laxed "asymmetric O-terminated" surface where every 
second surface O atom is removed and the others occupy 
the same sites as in the bulk structure. In this struc- 
ture, displacements of cations parallel to the surface are 
found substantially to lower the surf ace energy. A simila r 
study of the KNb0 3 (110) surface l)Heifets et aZll2000a|) . 
which is the surface of most experimental interest for 
this 1+/5+ perovskite, showed very strong relaxations 
extending deep below the surface, consistent with sug- 
gestions that this surface has a complicated chemistry. 

Next, we consider results on the electronic structure 
of BaTiC>3 films, particularly the surface states. In the 
FLAPW study of Cohen, the band gap of the ideal slab 
is is found to be reduced from the bulk. A primarily O 
p occupied surface state on the TiC>2 surface was identi- 
fied at M ijCoherl 1996), with a primarily Ti d surface 
state near the bottom of the conduction band. Analysis 
of the ferroelectric BaO-terminated slab showed that the 
macroscopic field resulted in a small charge transfer to 
the subsurface Ti d states from the O p and Ba p states at 
the other surface, ma king the surfaces meta llic. Further 
study of this effect in iKrcmar and FuL 120031 showed that 
for symmetric TiC>2 9-layer slab the ferroelectric distor- 
tion similarly shifts the top surface Ti states and bottom 
surface O state toward the bulk midgap as in Fig. re- 
sulting in a small charge transfer and a metallic character 
for the surfaces. 



b. PbTiO?, In contrast to the numerous papers on cal- 
culations on the surfaces of BaTiC>3, there are rela- 
tively few for the related material PbTiC>3. Regard- 
ing surface relaxations and energies, it was found in 
iMever and Vanderbiltll20oH that the two compounds are 
quite similar. In perpendicularly polarized films, it seems 
that both terminations give ferroelectric films if the de- 
polarization field is compensated, consistently with the 
stronger ferroelectric insta bility of PbTiO.3 and the mi- 
croscopic model analysis of lGhosez and Rabell2000l the 
latter not including the effects of surface relaxation. Be- 
cause of the larger spontaneous polarization of PbTiC>3, 
it is not possible fully to compensate the depolarization 
field using the dipole-layer technique. 

There arc important differences between A-site Ba and 
Pb, which are evident even for the bulk. While the po- 
larization in BaTiC>3 is dominated by the Ti displace- 
ments, Pb off-centering contributes substantially to the 
spontaneous polarization of PbTiOs; this can be linked 
to the much richer chemistry of Pb oxides compared to 
alkaline-earth oxides. One downside is that it is more 
challenging to construct accurate interatomic potentials 
for perovskites wi t h Pb t han with alkaline earth A cations 
<|Sepliarskv et aZl |2004|) . In the surface, the character- 
istic behavior of Pb leads to an antiferrodistortive sur- 
face reconstruction of the (001) PbO-t e rmina ted surface 
ijBungaro and Rabel iMunkholm et al I 12002ft . Specifi- 
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FIG. 24 (a) Paraelectric phase energy-band structure of nine- 
layer slab of BaTiOa . The surface states in (or near) the band 
gap are highligh ted, (b) As in (a), for the ferroelectric nine- 
layer slab. From lKrcmar and FuL l20oH . 

cally, first principles calculations ijBungaro and Rabel) 
show that the reconstruction in the subsurface TiC>2 layer 
occurs only for the PbO termination and not for TiC"2 
termination, and also that if the Pb in the surface layer 
is replaced by Ba, the reconstruction is suppressed. 



c. SrBi'2Ta20g A first-principles study of an isolated 
Bi02-terminated slab of SrBi2Ta20g (SBT) one lattice 
co nstant thick (com position SrBi 2 Ta 2 0n was reported 
in iTsai et all 1200.11 Spin-polarized calculations showed 
that such a film would be ferromagnetic as well as fer- 
roelectric. This intriguing possibility suggests further in- 
vestigation. 



d. SrTiC>3 and KTaC>3 In addition to the work on true 
ferroelectrics, there has been quite a lot of interest in 
first-principles studies of surfaces and heterostructures 
of incipient ferroelectrics, mainly SrTiOs and, to a lesser 
extent, KTaOs. As these have closely related properties 
that can illuminate issues in the ferroelectric perovskites, 
we include a description of a few representative results 
drawn from the very extensive literature on this subject. 
Many of the studies of BaTi0 3 discussed above in- 
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eluded analogous calculations for SrTiC>3. As already 
noted in the discussion of PbTi0 3 , the (001) surface re- 
laxations and energies of the nonpolar slab are very sim- 
ilar for all three materials. First principles surface re- 
laxation for the SrO surface is reported (in units of a 
= 3.86A) for surface Sr, surface O, subsurface Ti and 
subsurface O as -0.057, 0.001, 0.012 and 0.0, respectively 
(Padilla and Vanderbihill998|) . to be compared with the 
values -0.071, 0.012, 0.016, 0.009 (computed using a shell 
model with the experimental lattice constant 3.8969 A). 
Similarly, first principles surface relaxation for the TiC>2 
surface is reported for surface Ti, surface O, subsurface 
Sr and subsu rface O as -0.034, -0.016 , +0.02 5 and -0.005, 
respectively l|Padilla and Vanderbiltt IT998I) . to be com- 
pared with the values -0.030, -0.017, +0.035 and -0.021. 
The average energy of the two surfaces is found to be 
1.26 eV per surface unit cell. A detailed comparison of 
various Hartree-Fock and density-functional implemen- 
tations showed generally goo d agree ment for the sur- 
face relaxation ijHeifets el all l2001b|) . Inward surface 
dipoles due to relaxation are found for both termina- 
tions, with the TiC>2 ter mination to be small er in mag- 
nitude than for BaTi0 3 (|Heifets el aZJ . l2000b^ . The pos- 
sibility of an in-plane ferroelectric instability at the sur- 
face was examined and it wa s found to be quite weak 
ijPadilla and Vanderbihlll99S|) . In these studies, the an- 
tiferrodistortive instability exhibited by bulk SrTiC>3 at 
low temperatures was suppressed by the choice of a lxl 
in-plane unit cell. The surface electronic band structures 
show a behavior highly similar to that of BaTiC>3 de- 
scribed above. 

For SrTiC>3, there is considerable evidence for a wide 
variety of surface reconstructions of varying stoichiom- 
etry, depending on conditions such as temperature and 
oxygen partial pressure as well as the relative chemical 
potentials of TiC>2 and SrO. Candidate structures can 
be obtained by creating vacancies on the surface (for 
example, missing rows o f oxygen) and adding adatoms 
l)Kubo and Nozovell2003f) . More drastic rearrangements 
of the surface atoms have also been proposed , for example 
a (2x1) Ti 2 3 reconstruction (jCastelj 1200^) and a (2x1) 
double-lay er TiQ 2 reconst r uction with edge-sharing TiC>6 
octahedra ijErdman el al .1 120021) . As in the case of semi- 
conductor surface reconstructions, first-principles calcu- 
lations of total energies are an essential complement to 
experimental structural determination, and can also be 
used to predict scanning tunnell ing microscope image s 
for comparison with experiment l)Johnston el all [2004) . 
Even so, there are still many open questions about the 
atomic and electronic structure of SrTiC>3 surfaces un- 
der various conditions. The same applies to KTaOs; 
the structures and lattice dynamics of a variet y of (lxl) 
and (2 x1) surface structures were studied in l)Li el all 
12003am . 

As previously mentioned, the theoretical and experi- 
mental literature on SrTiOa is so extensive that it would 
require a review paper in its own right to cover it fully. 
Since, if anything, we expect the ferroelectric instability 



in related systems such as BaTiC>3 and PbTiC>3 to make 
the physics more, not less, complicated, this suggests that 
we have only scratched the surface in developing a com- 
plete understanding of the surfaces of perovskite ferro- 
electrics and the resulting effects on thin film properties. 



4. Studies of individual heterostructures 

Now we turn to the description of studies of systems 
with two or more material components. The main struc- 
tural issues are the rearrangements at the interface, the 
change in electrical and mechanical boundary conditions 
felt by each constituent layer, and how these changes 
modify the ferroelectric instability exhibited by the sys- 
tem as a whole. This geometry also allows the calculation 
of band offsets and/or Schottky barriers, crucial in prin- 
ciple to understanding the electronic behavior (though 
with the caveat that the measured Schottky barrier in 
real systems will be strongly influenced by effects such as 
oxidation of the electrodes that are not included in the 
highly idealized geometries studied theoretically). The 
current state of knowledge, derived from experimental 
measurements, is described in Section. IIII.B.3I 

The first combination we discuss is that of a ferro- 
electric thin film with metallic electrodes. Transition 
metal interfaces with n onpolar BaO-term inated layers of 
BaTiC>3 were studied in lRao el aZ.L 119971 specifically sys- 
tems of 3 and 7 atomic layers of BaTiC>3, with lattice 
constant set to the bulk value of 4.00 A, combined with 
top and bottom monolayers of Ta, W, Ir and Pt rep- 
resenting the electrodes. The preferred absorption site 
for the metal atoms was found to be above the O site, 
with calculated metal oxygen distances ranging from 2.05 
A for Ta to 2.11 A for Pt. The BaTi0 3 slabs were as- 
sumed to retain their ideal cubic structure. Analysis of 
the partial density of states of the heterostructure shows 
that the Pt and Ir Fermi energies lie in the gap of the 
BaTiC>3 layer at 0.94 eV and 0.64 eV, respectively, above 
the top of the valence band (this is, fortunately, smaller 
than the underestimated computed gap of 1.22 eV for 
the BaTiC>3 slab). Using the experimental gap of 3.13 
eV, a Schottky barrier height of 2.19 eV for Pt and 2.49 
eV for Ir is thus obtained. Experimentally, however, the 
Schottky barrier is known to be substantially lower for Ir 
than for Pt, illustrating the limitations mentioned in the 
prev ious paragraph. 

In lRobertson and Che'nlll999l first principles calcula- 
tions of the charge neutrality levels were combined with 
experimental values of the band gap, electronic dielectric 
constant e,^, the electron affinity, and the empirical pa- 
rameter S, described above in Section III.B 3. Values for 
SrTi0 3 and PZT were reported for Pt, Au, Ti, Al and 
the conduction and valence bands of Si. 

To explore how the electrodes affect the ferroelectric 
instability of the film, Junquera and Ghosez considered 
a supercell of 5 unit cells of metallic SrRu03 and 2- 
10 unit cells of BaTiC>3, with a SrO/Ti02 interface be- 
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tween the SrRuC>3 and the BaTiC>3. A SrTiC>3 substrate 
was treated implicitly by constraining the in-plane lat- 
tice constant of the supercell to that of bulk SrTi03. For 
each BaTi0 3 thickness, the system was relaxed assum- 
ing a nonpolar state for the BaTi03, and then the energy 
of the bulk-like tetragonal distortion was computed as a 
function of overall amplitude of the distortion. Above a 
critical thickness of six unit cells, this distortion lowered 
the energy, demonstrating the development of a ferro- 
electric instability. As will be discussed further in section 
IIV.B.5I this finite size effect can be largely understood by 
considering the imperfect screening in the metal layers. 

Considerable first-principles effort has been devoted 
to investigating various aspects of epitaxial ferroelectric 
thin films on Si. As perovskite oxides cannot be grown 
directly on Si, an approach devel oped by McKee and 
Walker l)McKee et all Il998l l200l|) is to include a AO 
buffer layer, which apparently also results in the forma- 
tion of a silicide interface phase. The constituent layers 
of this heterostructure should thus be considered to be 
Si/ASi 2 /AO/AB0 3 . The full system has not been sim- 
ulated directly, but first-principles approaches have been 
used to investigate individual interfaces. The importance 
of relaxations, the additional role of the buffer layer in 
changing the band offset, and the analysis of electronic 
structure within the LDA are illustrated by the following. 

In jMcKee et all 1200$ . first-principles results are pre- 
sented for the atomic arrangements and electronic struc- 
ture in the Si/ASi2/AO system, in conjunction with an 
experimental study. A strong correlation is found be- 
tween the valence band offset and the dipole associated 
with the A-0 bond linking the A atom in the silicide to 
the O atom in the oxide, shown in Fig. 123 It is thus 
seen that the structural rearrangements in the interface 
are a key determining factor in the band offset. 

A detailed examination of the interface between the 
perovskite oxide and the alkaline oxide buffer layer, 
specific ally BaO/BaTiO.3 and SrO/SrTi03, was carried 
out in ijjunauera et all l2003|) . A periodic 1x1x16 su- 
percell was chosen with stacking of (001) atomic layers: 
(AO) n -(AO-Ti02)m> with n — 6 and m = 5. Two mir- 
ror symmetry planes were fixed on the cental AO and 
BO2 layers, and the in-plane lattice constant chosen for 
perfect matching to the computed LDA lattice constant 
of Si (this epitaxial strain constraint is the only effect of 
the Si substrate included in the calculation). Relaxations 
within the highest-symmetry tetragonal space group con- 
sistent with this supercell were performed. Analysis of 
the partial density of states showed no interface-induced 
gap states. The main effect observed for relaxations was 
to control the size of the interface dipole, which in turn 
was found to control the band offsets, shown here in Fig- 
ure EH As in the studies described above, the conduc- 
tion band offset is obtained from the computed valence 
band offset by using the bulk experimental band gap. 
These results were combined with offsets reported for the 
other relevant interfaces to estimate band alignments for 
Si/SrO/SrTi0 3 /Pt and Si/BaO/BaTi0 3 /SrRu0 3 het- 
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FIG. 25 The image illustrates three layers of the alkaline 
earth oxide on the (001) face of silicon observed in a cross 
section at the [110] zone axis (blue, alkaline earth metal; yel- 
low, oxygen; and green, silicon). A distinct interface phase 
can be identified as a monolayer structure between the oxide 
and the silicon in which the charge density in interface states 
is strongly localized around the silicon atoms in the interface 
phase. The dipole in the ionic A-0 bond between the alka- 
line earth metal in the silicide and the oxygen in the oxide 
buffers the junction against the electrostatic polarization of 
the interface states localized on silicon. The electron den- 
sity of this valence surface state at the center of the Brillouin 
zone is shown with the purple isosurface (0.3 10 -3 e). From 
iMcKee et a/ll2003l 




FIG. 26 Schematic representation of the valence-band off- 
set (VBO) and the conduction-band offset (CBO) for 
BaO/BaTi0 3 (a) and SrO/SrTi0 3 (b) interfaces. E„, E c , 
and Egap* stand for the top of the valence band, the bot- 
tom of the conduction band, and the experimental band gap, 
respectively. Values for E„ , measured with respect to the 
average of the electrostatic potential in each material, are in- 
dicated. The solid curve represents the profile of the macro- 
scopic average of the total electrostatic potential across the 
interface. AV stands for the resulting lineup. The in-plane 
lattice constant was set up to the theoretical one of Si (5.389 ). 
The size of the superc ell corresponds to n=6 and m=5. From 
lJunauera et aUl2003l 
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erostructures, confirming that the AO layer introduces 
an electrostatic barrier of height greater than 1 eV. This 
is sufficient to eliminate the carrier injection from the 
Si into the conduction band states of the perovskite 
that would be occur if the two were in direct contact 
Robertson and CW ITflflft. 



5. First principles modelling: methods and lessons 

As discussed above, the analysis and prediction of the 
behavior of ferroelectric thin films and heterostructures 
can be carried out with direct first-principles simulations 
only for highly idealized configurations. However, it is 
possible to consider more complex and realistic situations 
by constructing models that incorporate certain physical 
ideas about the nature of these systems, with material- 
specific parameters determined by fitting to the results of 
first-principles calculations carefully selected for a com- 
bination of informativeness and tractability. This mod- 
elling approach also has the advantage of providing a 
conceptual framework for organizing the vast amount of 
microscopic information in large-scale first-principles cal- 
culations, and communicating those results, particularly 
to experimentalists. This will not diminish in importance 
even as such calculations become easier with continuing 
progress in algorithms and computer hardware. 

For successful modeling of measurable physical proper- 
ties, the film must be considered as part of a system (sub- 
strate+film+electrode) as all components of the system 
and their interaction contribute to determine properties 
such as the switchable polarization and the dielectric and 
piezoelectric response. We start by considering the class 
of first-principles models in which the constituent layers 
(film, superlattice layers, electrodes and substrate) are 
assumed to be subject to macroscopic electric fields and 
stresses resulting from the combination of applied fields 
and stresses and the effects of the other constituents, with 
the responses of the layers being given to lowest order by 
the bulk responses. For systems with constituent layer 
thicknesses as low as one bulk lattice constant, it seems 
at first unlikely that such an approach could be useful, 
but in practice it has been found to be surprisingly suc- 
cessful. 

One simple application of this approach has been 
used to predict and analyze the strain in nonpolar 
thin films and multilayers. In the construction of 
the reference structure for the AO/ABO3 interfaces in 
j.Tunauera et all |2003), macroscopic modelling of the 
structure with bulk elastic constants for the constituent 
layers yielded accurate estimates for the lattice constants 
along the normal direction. In cases of large lattice mis- 
match, very high strains can be obtained in very thin 
films and nonlinear contributions to the elastic energy 
can become important. These can be computed with 
a slightly more sophisticated though still very easy-to- 
implement method that has been developed to study the 
effects of epitaixial strain more generally on the structure 



and properties of a particular material, described next. 

As is discussed further in Section E] the effects of 
epitaxial strain in ultrathin films and heterostructures 
have been identified as a major factor in determining 
polarization-related properties, and have been the sub- 
ject of intense interest in both phenomenological and first 
principles modelling. In particular, for ferroelectric per- 
ovskite oxides it has long been known that there is a 
strong coupling between strain (e.g. pressure-induced) 
and the ferroelectric instability, as reflected by the fre- 
quency of the soft mode and the transition temperature. 
In both phenomenological and first-principles studies, it 
has become common to study the effects of epitaxial 
strain induced by the substrate by studying the struc- 
tural energetics of the strained bulk. Specifically, two 
of the lattice vectors of a bulk crystal are constrained 
to match the substrate and other structural degrees of 
freedom are allowed to relax, as described in the pre- 
vious paragraph. In most cases, these calculations are 
performed for zero macroscopic electric field, as would 
be the case for a film with perfe ct short-circuited elec- 
trodes. Indeed, it i s often the case d.Iunquera et qiJ . l2003t 
iNeaton and Rabel 120031: iPertsev et all Il998fl that the 
strain effect is considered to be the dominant effect of the 
substrate, which is otherwise not included (thus greatly 
simplifying the calculation). At zero temperature, the 
sequence of phases and phase boundaries can be readily 
identified as a function of in-plane strain directly through 
total-energy calculations of the relaxed structure subject 
to the appropriate constraints. Atomic-scale information 
can be obtained for the precise atomic positions, band- 
structure, and phonon frequencies and eigenvectors. The 
temperature axis in the phase diagram can be included 
by using effective Hamiltonian (or interatomic potential) 
simulations. Results for selected perovskite oxides are 
discussed in Section IV! a similar analysis w as reported 
for Ti02 in l(Montanari and Harrison! 12004^1 . 

This modelling is based on the assumption that the 
layer stays in a single-domain state. As discussed in 
Section EJ the possibility of strain relaxation through 
formation of multidomain structures must be allowed 
for. While this cannot be readily done directly in first- 
principles calculations, first-principles data on structural 
energetics for large misfit strains could be used to refine 
L andau parameters for use in calculations such as those 
inlAlpav and Rovtburdl ll998HBratkovsk v and LevanvkL 
l200li iLi et all l2003al ISpeck and PompeL 11994 The ef- 
fects of inhomogeneous strain due to misfit dislocations 
that provide elastic relaxation in thicker films have been 
also been argued to be significant. 

Next we consider the application of these "continuum" 
models to analyzing structures in which the macroscopic 
field is allowed to be nonzero. Macroscopic electrostatics 
is applied to the systems of interest by a coarse-graining 
over a lattice-constant-scale window to yield a value for 
the local macroscopic electric potential. Despite the fact 
that this is not strictly within the regime of validity of 
the classical theory of macroscopic electrostatics, which 
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requires slow variation over many lattice constants, this 
analysis turns out to be remarkably useful in the analysis 
of first principles results. In the simplest example, the po- 
larization of an polar BaTi03 slab (periodically repeated 
in a supercell with vacuum) is accurately reproduced us- 
ing bulk values for the bulk spontaneous polarization and 
electronic dielectr ic constant even for s labs as thin as two 
lattice constants IjJunauera and R abc). 

We have already mentioned in the previous section 
that perpendicular (to the surface) polarization can lead 
to a nonzero macroscopic field that opposes the po- 
larization (the depolarizing field). Unless compensated 
by fields from electrodes or applied fields, this strongly 
destabilizes the polarized state. In systems with two 
or more distinct constituent layers, this condition in 
the absence of free charge favors states with V • P = 
0. For example, in the first-principles calculations of 
short-period BaTi03/SrTiC>3 superlattices, the local po- 
larization along [001] is found to be quite uniform in 
the t wo layers ijjohnston and Rabel iNeaton and Rabel 
2003), though the in-plane component can be very dif- 
ferent l| Johnston and Rabel) . To the extent that layered 
ferroelectrics such as SBT can be treated in this macro- 
scopic framework, one similarly expects that polarization 
along c will tend to be energetically unfavorable, since the 
layers separating the polarized perovskite-like lay ers typ- 
ically have low polarizability ijFennie and Rabej) . This 
observation provides a theoretical framework for evalu- 
ating claims of large ferroelectric polar ization along c in 
layered compounds l)Chon et all l2002j) ; unless there is 
an unusually high polarizability for the non-perovskite 
layers, or a strong competing contribution to the energy 
due, for example, to the interfaces to help stabilize a high 
c polarization, other reasons for t he observations need to 
be considered llGarg et aZll2003ft . 

These considerations become particularly important 
for ultrathin films with metal electrodes. The limiting 
case of complete screening of the depolarizing field by 
perfect electrodes is never realized in real thin film sys- 
tems. The screening charge in real metal electrodes is 
spread over a characteristic screening length and is as- 
sociated with a voltage drop in the electrode. For thick 
films, this can be neglected, but the relative size of the 
voltage drop increases as the film thickness decreases. 
This has been identified as a dominant contribution to 
the relation between the applied field and the true fiel d 
in the film for the thinnest films ()Dawber et alll2003b|) . 
One way this shows up is in the thickness dependence of 
the apparent coercive field; it is found that the true coer- 
cive field scales uniformly down to the thinnest films. Ef- 
fects are also expected on the structure and polarization. 
While films with partial compensation of the depolariza- 
tion field may still exhibit a ferroelectric instability, the 
polarization and the energy gain relative to the nonpolar 
state are expected to decrease. Th is simple model was de- 
velop ed and successfully used in l)Junauera and Ghosed . 
2003) to describe the thickness dependence of the fer- 
roelectric instability in a BaTiC>3 film between SrRuC>3 



electrodes. This analysis identified the thickness depen- 
dence of the residual depolarization field as the prin- 
cipal source of thic kness dependence in this case. In 
iLichtensteieer et all the reduction of the uniform polar- 
ization by the residual field and its coupling to tetragonal 
strain was suggested to be the cause of the decrease in 
tetragonality with decreasing thickness of PbTiC>3 ultra- 
thin films. 

It is well known that 180 degree domain formation pro- 
vides an effective mechanism for compensating the depo- 
larization field, and is expected to be favored when the 
screening available from electrodes is poo r or nonexistent 
(for e xample, on an insulating substrate l)Streiffer et all 
2002) ). Instability to doma i n for mation is discussed 
in iBratkovskv and Levanvukl 12000' as the result of a 
nonzero residual depolarization field due to the presence 
of a passive layer. Similarly, a phase transition from a 
uniform polarized state to a 180 degree domain state with 
zero net p olarization is exp ected to occur with decreasing 
thickness l|Junauera et alb . 

Despite the usefulness of macroscopic models, it should 
not be forgotten that they are being applied far outside 
the regime of their formal validity (i.e length scales of 
many lattice constants) and that atomistic effects can 
be expected to play an important role, especially at the 
surfaces and interfaces. The structural energetics could 
be substantially altered by relaxations and reconstruc- 
tions (atomic rearrangements) at the surfaces and inter- 
faces. These relaxations and reconst ructions are also ex- 
pected to couple to the polar ization l|Bungaro and Rabel 
iMever and Vanderbiltl 12001^) . with the possibility of ei- 
ther enhancing or suppressing the switchable polariza- 
tion. The surfaces and the interfaces will also be primar- 
ily responsible for the asymmetry in energy between up 
and down directions for the polarization. For ultrathin 
films, the surface and interface energy can be important 
enough to dominate over elastic energy, leading to a pos- 
sible tradeoff between lattice matching and atomic-scale 
matching for favorable bonding at the interface. These 
surface and interface energies could even be large enough 
to stabilize non-bulk phases with potentially improved 
properties. This should be especially significant for in- 
terfaces between unlike materials. 

Electronic states associated with surfaces and inter- 
faces will also contribute to determining the equilibrium 
configuration of electric fields and polarizations. In the 
simple example of periodically repeated slabs separated 
by vacuum, as the slab gets thicker, a breakdown is ex- 
pected where the conduction band minimum on one sur- 
face of the slab falls below the valence band maximum on 
the other. In this case, charge will be transferred across 
the slab, with the equilibrium charge (for fixed atomic po- 
sitions) being determined by a combination of the macro- 
scopic electrostatic energy and the single particle density 
of states. This tends partially to screen the depolariza- 
tion field. The role of interface states in screening the 
depolarization field in th e film has be en discussed in a 
model for BaTiC>3 on Ge (R einer et alb . The presence of 
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surface and interface states can be established by exam- 
ination of the bandstructure and PDOS, as discussed in 
the previous section. 

Finally, we turn our attention to an analysis of what 
the discussion above tells us about finite-size effects in 
ferroelectric thin films. We have seen that many factors 
contribute to the thickness dependence of the ferroelec- 
tric instability: the thickness dependence of the depolar- 
ization field, the gradual relaxation of the in-plane lattice 
constant from full coherence with the substrate to its bulk 
value and the changing weight of the influence of surfaces 
and interfaces. The "true" finite-size effect, i.e the mod- 
ification of the collective ferroelectric instability due to 
the removal of material in the film relative to the infinite 
bulk, could possibly be distangled from the other factors 
by a carefully-designed first principles calculation, but 
this has not yet been done. We speculate that this effect 
does not universally act to suppress ferroelectricity, but 
cou ld, depending on the ma terial, enhance ferroelectric- 
ity l)Ghosez and Rabel l200fl . 



6. Challenges for first principles modelling 

First principles calculations have advanced tremen- 
dously in the last decade, to the point where systems 
of substantial chemical and structural complexity can be 
addressed, and a meaningful dialogue opened up between 
experimentalists and theorists. With these successes, the 
bar gets set ever higher, and the push is now to make the 
theory of ferroelectrics truly realistic. The highest long- 
term priorities include making finite temperature calcu- 
lations routine, proper treatment of the effects of defects 
and surfaces, and the description of structure and dy- 
namics on longer length and time scales. In addition, 
there are specific issues that have been raised that may 
be addressable in the shorter term through the interac- 
tion of theory and experiment, and the rest of this section 
will highlight some of these. 

Many applications depend on the stability of films with 
a uniform switchable polarization along the film nor- 
mal. This stability depends critically on compensation of 
the depolarization field. Understanding and controlling 
the compensation mechanism(s) are thus the subjects of 
intense current research interest. There are two main 
classes of mechanism: compensation by "free" charges (in 
electrodes/substrate or applied fields) and compensation 
by the formation of polarization domains. On insulat- 
ing substrates, this latter alternative has been observed 
and character i zed in ultrathin films l)Fong et all 120041 
IStreiffer et fl*ll2002ft : it has been proposed that domain 
formation will occur in films on conducting substrates at 
very low thicknesses as well where the finite screening 
length in realis tic electrodes in hibits that mechanism of 
compensation ijjunauera et alj} . The critical thickness 
for this instability depends on the domain wall energy. 
This is expected to be different in thin films than in bulk, 
one factor being that the bulk atomic plane shifts across 



the domain walls. 

Compensation of the depolarization field by free 
charges appears to be the dominant mechanism in films 
on conducting substrates (even relatively poor conduc- 
tors), with or without a top electrode. In the latter case 
there must be free charge on the surface; the challenge is 
to understand how the charge is stabilized. There are also 
unresolved questions about how the charge is distributed 
at the substrate-film interface, and how this couples to 
local atomic rearrangements. Asymmetry of the com- 
pensation mechanism may prove to be a significant con- 
tribution to the overall up-down asymmetry in the film 
discussed in the previous section. A better understanding 
could lead to the identification of system configurations 
with more complete compensation and thus an enhance- 
ment of stability. 

The study of the behavior of ferroelectrics in applied 
electric fields also promises progress in the relatively near 
future. Recently with the solution of long-standing ques- 
tions of principle, it has become possible to perform DFT 
calculations f or cry stalline solids in finite electric fields 
l)Souza et all 12002). In ferroelectrics, this allows the in- 
vestigation of nonlinearities in structure and polarization 
at fields relevant to experiments, and the possibility of 
more accurate modelling of constituent layers of thin film 
and superlattice systems subject to nonzero fields. It is 
also of interest to ask what the intrinsic breakdown field 
would be in the absence of defects, though the question 
is rather academic with respect to real systems. 

The nonzero conductivity of real ferroelectrics becomes 
particularly important for thinner films, both since a 
higher concentration of free carriers is expected to be 
associated with characteristic defects in the film, and 
because a given concentration of free carriers will have 
a more significant impact as thickness decreases. Free 
carriers can at least partially screen macroscopic elec- 
tric fields. At the macroscopic level, the concepts of 
band bending and space charge arising in semiconductor 
physics can be applied to thin film ferroelectrics, while 
a correct atomic-scale treatment of this effects could be 
important to describing the behavior of ultrathin films. 

The physics of switching presents a significant chal- 
lenge, requiring description of structure and dynamics 
on long length and time scales. The questions of what 
changes, if any, occur in switching as films become thin- 
ner continue to be debated. The possibility of switch- 
ing as a whole rather than via a domain- wall mechanism 
has b een raised for ultrathin films of PVDF ijBune et all 
1998), while a different interpretation has been offered in 
Dawber et all l2003bl Some progress has been made us- 
ing interatomic potentials for idealized defect-free films, 
though real systems certainly are affected by defects re- 
sponsible for such phenomena as imprint and fatigue. 
Ongoing comparison of characteristics such as coercive 
fields, time scales, material sensitivity, and thickness de- 
pendence of domain wall nucleation, formation energy, 
and motion with experimental studies promise that at 
least some of these issues will soon be better understood. 
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To conclude this section, we emphasize that it is not 
very realistic to expect first-principles calculations quan- 
titatively to predict all aspects of the behavior of chem- 
ically and structurally complex systems such as ferro- 
electric thin films, although successful predictions should 
continue to become increasingly possible and frequent. 
Rather, the quantitative microscopic information and the 
development of a useful conceptual framework contribute 
in a close interaction with experiment to build an under- 
standing of known phenomena and to propel the field 
into exciting new directions. 



V. STRAIN EFFECTS 

Macroscopic strain is a important factor in determin- 
ing the structure and behavior of very thin ferroelectric 
films. The primary origin of homogeneous film strain is 
lattice mismatch between the film and the substrate. In 
addition, defects characteristic of thin films can produce 
inhomogeneous strains that can affect the properties of 
thicker relaxed films of technological relevance. Because 
of the strong coupling of both homogeneous and inhomo- 
geneous strains to polarization, these strains have a sub- 
stantial impact on the structure, ferroelectric transition 
temperatures and related properties such as the dielectric 
and piezoelectric responses, which has been the subject 
of extensive experimental and theoretical investigation. 

The largest effects are expected in coherent epitaxial 
films. These films are sufficiently thin that the areal elas- 
tic energy density for straining the film to match the 
substrate at the interface is less than the energy cost 
for introducing misfit dislocations to relax the lattice pa- 
rameters back towards their unconstrained equilibrium 
values. (We note that for ultrathin films, the relaxed in- 
plane lattice constant will not in general not be the same 
as the bulk lattice constant, and the form er is m ore ap- 
propriate for computing lattice mismatch ijRabeT) .) Very 
high homogeneous strains, of the order of 2%, are achiev- 
able. For example, BTO films on STO, with a bulk mis- 
match of 2.2%, remain coh erent in equilibr ium up to a 
critical thickness of 2-4 nm ijSun et all l2004|) . With low- 
temperature growth techniques, the formation of misfit 
dislocations is kinetically inhibited and coherent films can 
be grown to thicknesses of t wo to three times the critical 
thickness llChoi et a/1 12004). Even these films, however, 
are much thinner than the minimum 120 nm thickness for 
films used in contemporary applications, and thus much 
of the discussion in this section is at present primarily of 
fundamental rather than technological interest. 

The structure of a coherent film can be a single-crystal 
monodomain structure, a polydomain structure, or even 
possibly multiphase. We discuss the simplest single- 
crystal monodomain case first. The phase diagram as 
a function of in-plane strain will in general include lower 
symmetry phases due to the symmetry-breaking charac- 
ter of the epitaxial constraint. A nomenclature for these 
phases of perovskites on a surface with square symmetry 




Misfit Strain u m , 10 -3 

FIG. 27 Phase diagram of a (001) single-domain PbTi03 thin 
films epitaxially grown on different cubic substrates providing 
various misfit strains um in the heterostructures. The second- 
and first-order phase tra nsitions are shown b y thin and thick 
lines, respectively. From lPertsev et al.lll99Sl 



(e. g. a perovskite (00 1) surface)) has been established 
in lPertsev et a/.lll998l For example, a ferroelectric per- 
ovskite rhombohedral phase will be lowered to monoclinic 
symmetry (called the r-phase). For highly compressive 
in-plane strains, coupling between strain and the polar- 
ization tends to favor the formation of a tetragonal phase 
with polarization along c (the c-phase) for highly com- 
pressive strains. Conversely, highly tensile strains lead to 
an orthorhombic phase with polarization along the cube 
face diagonal perpendicular to the normal (the aa-phase) 
or along the in-plane cartesian direction (the a phase). 
As a result of the added constraint, it is possible in prin- 
ciple to stabilize perovskite-derived phases not observed 
in bulk, for example, the monoclinic r-phase and the or- 
thorhombic aa-phase in PbTi03. This mechanism also 
plays a role in the more general pheno menon of epitaxial 
stabi lization, discussed, for example, in lGorbenko et al\ 
120021 

As discussed in Section TlV.B. 51 theoretical analysis of 
the in-plane strain phase diagrams has focused on iso- 
lating the effects of strain by computing bulk single- 
crystal monodomain phase diagrams under the epitaxial 
constraint and zero macroscopic electric field, using phe- 
nomenological Landau theory or first principles methods. 
Phenomenological analysis based on Landau-Devonshire 
theory for a number of perovskite o xides has been 
presented by Pertsey and co-work e rs ijK oipkjiar e t all 
l200lt iPertsev et aR I2003L l2000albl 119981 Il999j) . with 
temperature-strain diagrams for BaTiOs, PbTiC>3 (see 
Fig. EJ, SrTi0 3 , and Pb(Zr,Ti)0 3 (PZT), the latter 
generalized to include nonzero stress. First-principles 
methods have been used to construct a t e mpera ture- 
strain diagram for BaTiC>3 l)Dieguez et all l2004[) and 
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FIG. 28 Polarization of SrTi03 as a function of in-plane 
strain. Solid circles and squares denote polariz ation along 
[001] and [110], respectively. From lAntons et all 



a zero-temper ature strain diagram for PbTi03 and 
ordered PZT ijBungaro and Rabel. l2004|) and SrTi0 3 
l)Antons et alb . These theoretical phase diagrams have 
some notable features. In particular, compressive in- 
plane strain is found to elevate the ferroelectric (c)- 
paraclcctric transition temperatures in BaTiOa and 
PbTiC>3, and tensile in-plane strain elevates the ferroelec- 
tric (aa)-paraelectric transition temperatures. In both 
cases, the transition is second order, in contrast to the 
first-order transition in bulk. To eliminate a possible 
source of conf usion, we comment t hat zero misfit strain 
as defined in (jPertsev et all Il998j) is not equivalent to 
an unconstrained film (the low-temperature bulk phases 
are in general not cubic, and the constrained dielectric 
and piezoelectric responses are clamped, as will be dis- 
cussed further below). The nearly vertical morphotropic 
phase boundar y characteristic of bulk PZT is subs tan- 
tially modified l)Li et all l2003bt iPertsev et a/.l.l2003|) . In 
SrTiC>3, ferroelectricity is found to be induced by both 
sufficiently compressive and tensile strains, with a corre- 
sponding direction for the spontaneous polarization (c - 
type and aa-type) <|Antons et all IPertsev et aILl2000a(l . 
as shown in Fig. |^H1 The enhancement of the polariza- 
tion in the c-phase by c ompressive in-plane strain has 
been noted for BaTiOa IINeaton and Rabei . 120031 ) and 



PZT jPertsev et all l2003fT For both strained SrTi0 3 
and strained BaTiC-3, the ferroelectric T c 's are predicted 
to increase as the strain magnitudes increase (Choi 2004, 
Haeni 2004), as shown in Fig. HO 

The use of phenomenological bulk Landau parame- 
ters yields a very accurate description for small strains 
near the bulk T c . However, different parameter sets 
have been shown, for example in the case of BaTiC>3, 
to extrapolate to qu alitatively different phase d iagrams 
at low temperatures ijPertsev et all Il99a. Il999j) . Quan- 
titatively, the uncertainty in predicted phase boundaries, 
produced by the fitting of the Landau theory parameter 
inc reases with increa sing misfit strain as shown in Fig. 
1551 ( Choi etaR 120041) . The Landau analysis is thus well- 
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FIG. 29 Expected T c of (a) (001) BaTiQ 3 (fromlChoi et all 
(2004)) and (b) (001) SrTiOs (from lHaeni et alH2004) ) based 
on thermodynamic analysis. The range of transition repre- 
sents the uncertainty in the predicted T c resulting from the 
spread in reported property coefficients. 



complemented by first-principles calculations, which can 
provide very accurate results in the limit of zero temper- 
ature. In the case of BaTiC^ , the ambiguity in the low- 
temperature phase diagram (IPertsev et all I1998L Il999|) 
has been resolved in this way ipieeuez et all |2004|) ~ m 
the cas e of PbTiOa, the phenome nological result is con- 
firmed l|Bungaro and Rabel 120041) . 

The epitaxial-strain induced changes in structure and 
polarization are also expected to have a substantial ef- 
fect on the dielectric and piezoelectric responses. While 
overall, the dielectric and piezoelectric re sponses should 
be re duced by clamping to the substrate l|Canedv et all 
1200(1 iLi et flfl . l200ll) . these response s will tend to diverge 
near s econd-order phase bo undaries l|Bungaro and Rabel 
12004 IPertsev et all |2003|) Responses will also be large 
in phases, such as the r-phase, in which the direction 
of the polarization is not fixed by symmetry, so that 
an applied field or stres s can rotate the polarization 
l)Wu and Krakauerl Il992|) . This polarization rotation 
has been identified as a key mechanism in the colossal 
piezoresponse of single-crystal relaxors l|Fu and Cohenl 
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FIG. 30 Inverses of the eigenvalues of the phonon contribu- 
tion to the dielectric tensor as a function of the in-plane lat- 
tice constant for the [001]- (PbTiQ3) 1 (Pb ZrQ3) 1 superlattice. 
The lines are a guide to the eye. From IfTungaro and Rabel 
12004 



l2000t iPark and Shroud. Il997|) . The sensitivity of the 
zero-field responses should also be reflected in the non- 
linear response; thus the electric-field tunability of the 
dielectric r esponse can be ad justed by changing the mis- 
fit strains l|Chen et a/.Ll2003[) . 

While the phase diagrams thus derived are quite rich, 
even the optimal single-crystal monodomain structure 
may b e unfavorable with respect to formation of polydo- 
main l|Bj^atkoyj$k y_jm^ l200lt iRovtburd et all . 

120011 ISpeck and PompeLll994|) or multiphase structures, 
which allow strain relaxation on average and reduce elas- 
tic energy. The evaluation of the energies of polydomain 
structures requires taking both strain and depolarization 
fields into account. A recent d i scussio n of PbTiC>3 using 
a phase-field analysis l|Li et adl2003b|) suggested that in- 
cluding the possibility of polydomain structure formation 
significantly affects the phase diagram for experimentally 
relevant misfit strains and temperatures, as shown in Fig. 
|3"T1 For very thin films, additional effects associated with 
the interface between the film and substrate are also ex- 
pected to contribute significantly. For example, although 
polydomain formation can accommodate misfit strain av- 
eraged over different variants, each domain will be mis- 
matched to the substrate at the atomic level, with a cor- 
responding increase in interface energy. Furthermore, the 
energy of a domain wall perpendicular to the substrate 
will be higher than the energy of the corresponding wall 
in the bulk, due to the geometrical constraint on the al- 
lowed shifts of the atomic planes across the doma i n wal l 
(as found for the bulk in iMever and Vanderbifd . l2002f) 
imposed by the planar interface. Different domain walls 
will in general be affected differently by the constraint, 
possibly changing the relative energy of different polydo- 
main configurations. 

With recent advances in thin film synthesis, it is pos- 
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FIG. 31 (a)Phase diagram of PZT film under in-plane ten- 
sile strain of 0.005 obtained using thermodynamic calcula- 
tions assuming a single - domain state. There are only 
two stable ferroelectric phases. The solid lines represent 
the boundaries separating the stability fields of the paraelec- 
tric and ferroelectric phases, or the ferroelectric orthorhom- 
bic and distorted rhombohedral phases, (b) Superposition 
of the phase diagram of a PZT film under in-plane tensile 
strain of 0.005 from the phase-field approach (scattered sym- 
bols) and from thermodynamic calculations assuming a sin- 
gle domain (solid lines). There are three stable ferroelectric 
phases: tetragonal-" square," orthorhombic-"circle," and dis- 
torted rhombohedral-" triangle" according to the phase- field 
simulations. The scattered symbols represent the ferroelectric 
domain state obtained at the end of a phase-field simulation. 
The shaded portion surrounded by the scattered symbols la- 
bel the stability regions of a single ferroelectric phase, and the 
nonshaded re gion shows a mix ture of two or three ferroelectric 
phases. From lLi et aZ.ll2003bl 



sible to grow and characterize high-quality films that 
are sufficiently thin to be coherent or partially relaxed. 
Here, we give a few examples of experimental observa- 
tions of changes in structure and polarization in very 
thin films. T c elevation in s trained PbTiO s films has 
been reported and an alyzed in lRossetti et all l)l99l|) and 
IStreiffer et all l)2002|) . The strain- induced r-phase has 
been observed in PZ T films thinner than 150 nm on Ir- 
electroded Si wafers ijKelman et all . 12002^1 . An antifer- 
roelectric to ferroelectric transition has been observed 
in thin films of PbZr0 3 l)Avvub et all . Il998f) . though 
whether it is induced by strain or some other thin-film 
related effect is considered an open question. Large po- 
larization enhancements have been observed in epitax- 
ially strained BaTi0 3 films llChoi et all . 12004^ ■ Most 
dramatically, room-temperature ferroelectricity has been 
achieved for SrTi0 3 un der biaxial tensile s train induced 
by a DySc0 3 substrate ijHaeni et all . l2004|) . For thicker 
films, observation of polydomain structures is reported in 
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FIG. 32 Evolution of a x as a function of film 
thickness for Bao.eSro^TiOa thin films grown on 
0.29(LaAlO 3 ):0.35(Sr 2 TaAlO 6 ) substrates. .Also shown 
is the theoretical curve, given by the open circles. The 
straight dashed line represents the latt ice parameter of the 
ceramic target (a = 0.39505 nm). F rom ICanedv et all I2000 . 



iRovtburd et all l200ll It has been observed that domain 
formation may be suppressed by ra pid cooling through 
the transition l|Ramesh et a/.lll9 93). 

For thicker films grown at high temperature, misfit 
dislocations will form at the growth temperature par- 
tially or completely to relax misfit strain. The degree 
of relaxation increases with increasing thickness, un- 
til, for thick enough films, the epitaxial strain is neg- 
ligible. This behavior has bee n studied theoretically 
(|Matthews and Blakesleel Il974l) and observed experi- 
mentally, for example as in Fig. [321 Additional strain 
can arise during cooling from the growth temperature if 
there is differential thermal expansion between the film 
and the substrate and the formation of misfit dislocations 
is kinetically inhibited. A detailed theoretical study of 
strain relaxation in epitaxial ferroelectric films, with dis- 
cussion of the interplay of misfit dislocations, mixed do- 
main formation and de polarizing energy, was un dertaken 
by Speck and Pompe ijSpeck and Pompel 11994) - It was 
assumed that for rapid cooling from the growth tempera- 
ture, the effect of misfit dislocations can be incorporated 
by using an effective substrate lattice parameter, while in 
the limit of slow cooling, the system optimally accommo- 
dates misfit strain with dislocations. (This assumption 
is valid for films with thickness of order 1 /i, while the 
treatment needs to be slightly modified for intermediate 
thicknesses where the equilibrium concentration of misfit 
dislocations leads to only partial strain relaxation) . Elas- 
tic domains form to relax any residual strain. Below T c , 
depolarizing energy can change the relative energetics of 
different arrangements of polarized domains and misfit 
dislocations. It was suggested that the electrostatic en- 
ergy should be more of a factor for smaller tetragonality 
systems (BaTi0 3 vs PbTi0 3 ) where the strain energy is 
less, though this could at least be partially balanced by 
the fact that the polarization is smaller as well. A typical 
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FIG. 33 Coherent temperature dependent domain stability 
map for PbTiOa including the cubic lattice parameter for 
several common single crystal oxide substrates. The misfit 
strains for epitaxial growth of PbTi03 at 600° C are included 
in the insets along with t he critical thickness h for misfit dis- 
location formation. From ISpeck and Pompelll99l 



coherent diagram is shown in Fig. 

With transmission electron microscopy (TEM) it is 
possible to make detailed studies of the types and ar- 
rangements of misfit dislocations in perovs kite thin films . 
Rece nt st udies of high-qua lity films include fSuzuki et all 
EH and [S un et al 1 12004 While strain-relaxing defects, 
such as misfit dislocations, reduce or eliminate the elas- 
tic energy associated with homogeneous strain, these and 
other defects prevalent in films do generate inhomoge- 
neous strains. As mentioned at the beginning of this sec- 
tion, the inhomogeneous strains couple strongly to the 
polarization, and it has been shown by phenomenologi- 
cal analysis ijBalzar and Pooal 12004 iBalzar et all [2002) 
that their effects on T c can be significant. They have 
also been argued to contribute to the degradation of the 
dielectric response in thin films relative to bulk values 
ICanedv et «iJ . l2000D . 

Strains and their coupling to polarization are also cen- 
tral to the properties exhibited by short-period superlat- 
tices of lattice-mismatched constituents. As the result 
of recent work on artificial supcrlattices of ferroelectric 
materials, there are some indications that improved ferro- 
electric properties and/or very large dielectric constants 
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BaTi0 3 /SrTi0 3 superlattices 
lattice-matched to a SrTi0 3 substrate, the compressive 
in-plane strain on the BaTi0 3 layer substantially raises 
its polarization. Theoretical studies suggest that the 
SrTi0 3 layer is polarized (and the polarization in the 
BaTi0 3 layer is reduced) by electrostatic energy consid- 
erations, which favor continuity of the component of the 
polarization along the normal. Overall the polarization is 
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enhanced above that of bulk BaTi03, though not as high 
as that of a pure coherent BaTi03 film if it were possi- 
ble to suppress the formation of strain-relaxing defects. 
While the natural lattice constant of BaTiOa/SrTiOa is 
intermediate between the two endpoints, so that on a 
SrTiC>3 substrate the superlattice is under compressive 
in-plane stress, it has been suggested that the multilayer 
structure tends to inhibit the formation of misfit disloca- 
tions so that a thicker layer of coherent superlattice ma- 
terial can be grown. As the superlattice material thick- 
ness increases, there will be strain relaxation via misfit 
dislocations and the in-plane lattice constant should in- 
crease, putting the SrTiOs layer under in-plane tensile 
strain. In this case the SrTi03 layer is observed to have a 
component of pol arization along [110] ij.Tiang et a/.ll2003i 
iRios et ah . l2003|) , consiste nt with theoretical studies of 
epitax ially strained SrTiC-3 l)Antons et a^llPertsev et ah . 
l2000a|) and of the BaTiOa/SrTi O^ superlattice with e x- 
panded in-plane lattice constant II Johnston and Ra be). 

The real appeal of short-periodicity ferroelectric 
multilayers is the potential to make "new" artificially 
structured materials with properties that could open the 
door to substantial improvements in device performance 
or even radically new types of devices. Perovskites 
are particularly promising, as individual materials 
possess a wide variety of structural, magnetic, and 
electronic properties, while their common structure 
allows matching at the interface to grow superlattices. 
Beyond the prototypical example of BaTiC>3/SrTi03 
discussed in the previous paragraph, there has been 
work on other combinations such as KNbOs/KTaOs 
llChristen et all Il996t ISeoliarskv eTah 120021 12001 ; 
Sigman et all 120021) . Pb TiOs /SrTiOs ll.Iiang et al . 
1999}) . PbTi(WPbZrOs l|Bungaro and RabeL I200S , 
2004) . 
1999), 



VI. NANOSCALE FERROELECTRICS 



A. Quantum confinement energies 



Lan.fiSrn^MnOs/ Lan.fiSrn^F cOatlzumi et ... .. 
CaMn Qcj /CaRuOa |Takahashi et all *2001), 
LaCrO^- LaFeO^Ueda 6^11119981 119994 .and LaFe0 3 - 
LaMnCh l)Ueda et all Il999b() . In nearly all cases, 
strain plays an important role in understanding the 
aggregate properties of these short-period multilayers 
and superlattices. In addition to lattice mismatch, 
the layers also interact through the mismatch in 
polarization along the layer normal, which leads to 
mutual influences governed by considerations of elec- 
trostatic energy and nonzero macroscopic electric 
fields. With three or more constituents, it is possible 
to break inversion symmetry to obtain superlattice 
materials with possibly favorable piezoelectric prop- 
erties. This i dea was first proposed theoretically 
(|Sai et all l2000j) . leading to experimental stu d ies o f 
CaTi0 3 /SrTi0 3 /BaTi03 l|Warusawithana et all 120031) 
and LaA103/(La , Sr)Mn03/SrTiCh dKimoto et all 
12004 iQgawa et all 120031 lYamada et oil l2002|) . Also, 
perovskite superlattices combining ferroelectric and 
ferromagnetic layers offer a path to the development of 
multiferroic materials. The identification, synthesis and 
characterization of further combinations remains the 
subject of active research interest. 



Confinement energies are a trendy topi c in nanoscale 
semiconductor microelectronics devices. ijPetroff et all 
1200 1|) The basic idea is that in a system in which the elec- 
tron mean free path is long with respect to the lateral di- 
mension^) of the device, a quantum-mechanical increase 
in energy (and of the bandgap) in the semiconductor will 
occur. In general confinement energies exist only in the 
ballistic regime of conduction electrons, that is, where the 
electron mean free path exceeds the dimensions of the 
crystal. This usually requires a high-mobility semicon- 
ductor at ultra-low temperatures. Such effects are both 
interesting and important in conventional semiconduc- 
tors such as Si or Ge, GaAs and other III-Vs, and perhaps 
in II- Vis. However, despite the fact that the commonly 
used oxide ferroelectrics are wide-ba ndgap p-type semi- 
cond uctors (3.0 eV < E s < 4.5 eV). (|Waser and Smvthl 
1996) neither their electron nor hole mean free paths are 
sufficiently long for any confinement energies to be mea- 
sured. Typically the electron mean free p ath in an ABO3 
ferroelectric perovskite is 0.1 to 1.0 nm. ljDekker I Il9&4|) 
depending on applied electric field E, whereas the de- 
vice size d is at least 20 nm. Therefore any confinement 
energy (which scales as d~ 2 ) might be a meV or two, 
virtually unmea s urable, despite a few published claims 
l|Yu et a/.til997MKohiki et adl2000() . l|J.F Scottl l2000|) 
reporting extraordinarily large effects. In the case of 
Bi 2 03 and SrBi 2 Ta20g (SBT) these effects ma y arise 
from t wo-phase regi o ns at the sample surfaces. l|Zhou I 
ll992|) . (|Switzer et all Il999|) This is theoretically inter- 
esting and very important from an engineering device 
point of view; if it were not true the contact poten- 
tial at the electrode interface in a 1T-1C device, or at 
the ferroelectric-Si interface in a ferroelectric-gate FET, 
would depend critically on the cell size, which would add 
a very undesirable complication to device design. 



B. Coercive fields in nanodevices 



One of the most pleasant surprises in the research 
on small-area ferroelectrics is the observation, shown in 
Fig. EH , t hat the coerc ive field is independent of lateral 
area. ljAlexe et aZJ . Il999l) Coercive fields in nanophase fer- 
roelectric cells have gener ally been measured via a tomic 
force microscopy f AFM) . ijGruverman et all 1996) Do- 
main structures, polarization and coercive fields of 
nanoscale particles of BaTiOs hav e been studied the oret- 
ically using interatomic potentials ( Stachiotti , 20 Q4|) and 
a firs t-principles effective Hamiltonian l)Fu ImdBeUaichel 
120031) 
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FIG. 34 Lack of significant depende nce of coerci v e field on 
lateral area in nanoscale ferroelectrics \ Alexe et, all \199i) 



C. Self Patterned nanoscale ferroelectrics 

One approach to producing nanoscale ferroelectrics is 
to attempt to produce self patterned arrays of nanocrys- 
tals, in which ordering is produced by interactions be- 
tween islands through the substrate. This approach could 
be used to produce arrays of metallic nanoelectrodes on 
top of a ferroelectric film or alternatively arrays of crys- 
tals from the ferroelectric materials th emselves. The first 
scheme was suggested by Alexe et al. (jAlexe et aZJ . Il99^ 
who found that a bismuth oxide wetting layer on top of 
a bismuth titanate film formed an array of metallic bis- 
muth oxide nanocrystals on top of the film, which were 
partially registered along the crystallographic directions 
of the underlying substrate (Fig. I35fl . These nanocrys- 
tals were u sed successfu l ly as e lectrodes to switch regions 
of the filni l|Alexe et all Il999j) . In the second approach 



one might use a material such as PbTiC>3 on a SrTiC>3 
substrate, which was first demonstrated to form islands 
when grown e pitaxially at very t hin film thicknesses by 
Seifcrt et al. ijSeifert et all Il996|) In the context of self 
patterning of oxide materials a recent work by Vasco et 
al studies t he growth of self o rganised SrRu03 crystals 
on LaA10 3 l)Vasco et a/.ll2003[) . 

When small amounts of materials are deposited on 
substrates where there is some degree of mismatch be- 
tween the two materials, islands form and the repulsive 
interactions between them are mediated via st rain fields 
in th e substrate as first suggested by Andreev. l)Andreevl 
This idea has been develo ped into a detailed the- 
ory b y Shchukin and Bimberg;( Shchuki n and BimbereL 
1999) however this theory is a zero-temperature the- 
ory, whereas a thermodynamic theory is required to de- 
scribe the crystallization processes which occur at quite 
high temperatures. An extension of the theory to finite 
temperatures has been carried out by Williams and co- 
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FIG. 35 Sample of Alexe et al lAlexe et &\1\199A) . (a) TEM 

cross-section showing underlying layers and bismuth oxide na- 
noeelectrodes, (b) Semi-registered array of nanoelectrodes tak- 
ing their orientation from the underlying Si substrate 



workers . ijWilliams et al 1 l2000HRudd et all |2003) The 
chief result of this theory are the prediction of three dif- 
ferent kinds of structures (pyramids, domes and super- 
domes), a volume distribution for a particular species of 
structure, and a shape map to describe relative popula- 
tions of structures as a function of coverage and crystal- 
lization temperature. One interesting result from exper- 
iment is that similar shaped structures are observed in 
both the Volmer- Weber (VW) and Stranski-Krastanow 
(SK) growth modes, but on different size scales. In the 
work of Williams the thickness above which dome pop- 
ulations occur is of the order of 4-5 monolayers, cor- 
responding to the critical thickness for misfit disloca- 
tion s for Ge on Si(lOO). O n the other hand Capellini et 
al. ijCapellini et all Il997) studied via atomic force mi- 
croscopy the growth of Ge on Si(100) in the SK growth 
mode and found a much larger critical structure height 
of 50 nm at which dislocations were introduced and the 
structures changed from being pyramidal in geometry to 
domelike. The large increase in critical thickness is due 
to a substantial part of the misfit strain being taken up 
by the substrate in the S K growth mode, as descr i bed b y 
Eaglesham and Cerrulo (Eaglcs ham and Cerullol |l990|) . 
The description of self patterned ferroelectric nanocrys- 
tals by the models of Schukin and William s has recently 
been undertaken bv lDawber et al\ l)2003c^ . 

Prior to this two groups have grown PbTiOa nanocrys- 
tals on Pt/Si(lll) substrates to mea s ure size effects 
in fer roelectricity. ijRoelofs et all 12003(1 1 l)Shimizu et all 
120031) These works both show a lack of piezoresponsc 
in structures below 20nm in lateral size (Fig. I36fl . 
though we expect that this is connected to mechanical 
constraints rather tha n any fundamen tal limiting size for 
ferroelectric systems. IChu et al\ (|2004f) have highlighted 
the role that misfit dislocations can play in hampering 
ferroele ctricity in sma l l struc tures . Interestingly in the 
work of iRoelofs eTal\ i|2003ft and IShimizu et al\ l|2003ft 
because of the (111) orientation of their substrates, in- 
stead of square-based pyramids they obtain triangular 
based structures that display hexagonal rather than cu- 
bic registration (an analo gous result is observe d when 
Ge is grown on Si(lll). l|Capellini etail Il999fl ). The 
growth and analysis of PZT nanocrystals on SrTiOs has 
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(a) Topograpic image (b) Piezoresponse 
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FIG. 36 Roelof s~ et alJ \200A) (a) topographic image of grains 
from 100 nm down to 20 nm in lateral size (b) piezoresponse 
image of same grains showing the absence of peizoresponse for 
the grains below 20nm 

been carried out bv ISzafraniak et all l|2003l) 1. A review 
on size effects i n ferroelectric nanocry stals is currently in 
preparation bv lRuedieer et al\ l|2004|) . 

Although there is potential to produce self patterned 
arrays with greater registration by better choice of ma- 
terials and processing conditions our general conclusion 
is that highly registered memory arrays will not occur 
spontaneously in the absence of a pre-patterned field. 

D. Non-planar geometries:Ferroelectric nanotubes 

Almost all recent work on ferroelectric oxide films have 
involved planar geometries. However, from both a device 
engineering point of view and from theoretical consider- 
ations, it is now appropriate to analyze carefully non- 
planar geometries, especially nanotubes. 

Nanotubes made of oxide insulators have a variety 
of applications for pyroelectric detectors, piezoelectric 
ink-jet printers, and memory ca pacitors that cannot 
be fi l led by other nanotu bes llHerzog and Kattnerl 
1 19851) . ([Sakamaki et all 120011) JSaieey and Buschl 

l2002|) . jGnade et all l200(i|) . l|Averdung et all l200l[) . In 
the drive for increased storage density in FRAM and 
DRAM devices, complicated stacking geometries, 3D 
structures and trenches with high aspect ratios are 
also being investigated to increase the dielectric surface 
area. The integration of ferroelectric nanotubes into Si 
substrates is particularly important in construction of 
3D memory devices beyond the present stacking and 
trenching designs, which according to the international 
ULSI schedule 3 must be achieved by 2008. Template 
synthesis of nanotubes and wires is a versatile and in- 
expensive technique for producing nanostructurcs. The 



3 International Technology Roadmap for Semi- 
conductors (ITRS) 2002 (available at 
http://public.itrs.net/Files/2002Update/Home.pdf) 



size, shape and structural properties of the assembly are 
simply controlled by the template used. Using carbon 
nanotubes as templates, tubular forms of a number of 
oxides includin g V2OR, SiOo; AI2 O3 and Z1O2 have 
been generated ijPatzke et alJ.l2002|) . Much larger (>20- 
micron diameter) ferroelectric micro-tubes have been 
made by sputter depositio n around polyester fibres ([FoxL 
119951) . J Pokropivnvt l200l|) - Fox has made them from 
ZnO and PZT, with 23 /Ltm inside diameter, about lOOOx 
larger than the smallest nano-tubes reported in the 
present paper. Porous sacrificial templates as opposed 
to fibres have also been used. Porous anodic alumina 
has a polycrystalline structure with ordered domains of 
diameter 1-3 /im, containing self-organised 2D hexagonal 
tubul ar pore a rrays with an interpore distance of 50-420 
nm l)Li et all Il998|) . This nano-channel material can 
therefore be used as a template for individual nanotubes 
but is not suitable for making an ordered array of tubes 
over length scales greater than a few mm. Many oxide 
nanotubes, such as TiC>2, In 2 03, Ga 2 03, BaTiC>3 and 
PbTiC>3, as well as nanorods of M11O2, C03O4 and 
Ti02, have bee n made using porous alu mina membranes 
as tem plates (|Patzke et all l2002|) ]. IHernandez et al\ 
( 2002), used a sol-gel template synthesis route to prepare 
BaTiG-3 and PbTiC>3 nanotube bundles by dipping alu- 
mina membranes with 200 nm pores into the appropriate 
sol. The BaTiC>3 and PbTiO*3 nanotubes were shown 
to be cubic (paraelectric) and tetragonal (ferroelectric) 
by x-ray diffraction, although Raman studies indicated 
some non-centrosymmetric phase on a local scale in 
the BaTiC>3. Porous silic on materials a r e also available 
as suitable templates. iMishina et al\ (|2002l) used a 
sol-gel dipping technique to fill nanoporous silicon with 
a PbZri-zTizOs (PZT) sol producing nanograins and 
nanorods 10-20 nm in diameter. The presence of the 
ferroelectric PZT phase was shown by second harmonic 
generation (SHG) measurements. In this instance 
the p orous silicon does not ha ve a periodic array of 
pores llSmith and CollinsL |l992) and as in the case for 
those produced by Hernandez et al, we emphasise that 
those nanotubes are not ordered arrays, but instead 
spaghetti-like tangles of nano-tubes that cannot be 
used for the Si device embodiments. A second type of 
porous Si templates, however, consist of a very regular 
periodic array of pores with very high aspect ratios. By 
a combination of photolithography and electrochemical 
etching hexagonal or orthogonal arrays of pores with 
diameters 400 nm to a few mm and up to 100 ixm deep 
can b e formed in single cry stal Si wafers l|Schilling et all 
l200l|) . (|Ottow et aZ.l Il996^ . These crystals were origi- 
nally developed for application as 2D photonic crystals, 
but also find applications as substrates for templated 
growth and inte gration of oxides nanostructures with 
Si technology. iLuo et all l|2003f) recently used such 
crystals to produce individual, free standing PZT and 
BaTiC>3 ferr oelectric nanotu b es by a polymeric wetting 
technique. iMorrison et al\ l)2003j) described the use 
of liquid source misted chemical deposition (LSMCD) 
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to fill such photonic Si crystals with SBT precursor. 
During deposition, the SBT precursor was shown to coat 
the inside of the pores. After etching of the photonic 
crystal with pore diameter 2 /im for 30 seconds with aq. 
HF/HNO3 the interface between the Si substrate and 
SBT coating is dissolved, exposing the uniform SBT 
tube, Fig. 137b ,. The tube walls are very uniform with a 
thickness of ca. 200 nm. The same sample is shown in 
cross-sectional view after complete removal of the host 
Si walls between pores, Fig. ET7b . The result is a regular 
array of tubes attached to the host Si matrix only at the 
tube base. Although these tubes have suffered damage 
during handling, it is clear that the pores have been 
filled uniformly to the bottom, a depth of ca 100 /im. 
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FIG. 37 SEM micrograph indicating a plan view of a regular- 
array of SBT tubes in host silicon substrate with diameter 
ca. 2 fim and wall thickness ca. 200 nm (a). SBT tubes in 
cross sectional view indicating coating to bottom of pore (b). 
Micrograph of free-standing array of tubes with diameter ca. 
800 nm (c) and wall thickness < 100 nm (d). 

The second photonic crystal with pore diameter 800 
nm underwent fewer depositions and after etching re- 
vealed a regular array of uniform tubes of diameter 800 
nm, Fig. 137b . The wall thickness is uniform and < 100 
nm, Fig. I3TI 1. The tubes are ca. 100 fim long, completely 
discrete and are still attached to the host Si matrix, creat- 
ing a perfectly registered hexagonal array. Free-standing 
tubes may be produced by completely dissolving the host 
Si matrix. As yet, no one h as applied cylindrical e lec- 
trodes to the tubes; however. ISteinhart et all l|2003ft re- 
cently used porous anodic alumina templates to grow 
palladium nanotubes. Using a similar method it may 
be possible to alternately deposit Pd or Pt and SBT to 
produce a concentric electrode/FE/elctrode structure in 
each nanotube. The use of the photonic crystal template 
with a regular array of pores has significant benefits over 
other porous substrates in that the coatings/tubes pro- 
duced are also in a registered array ordered over several 
mm's or even em's. This facilitates addressing of such an 
array for device applications. DRAMs utilise high sur- 



face area dielectrics, and high aspect ratio SBT coatings 
such as these embedded in Si could increase storage den- 
sity. Current state-of-the-art deep trenched capacitors 
are 0.1 mm diameter by 6 mm deep, aspect ratio 60:1. 
Using SBT (or other FE oxide) nanotubes of wall thick- 
ness <100 nm, a trench (or array of trenches) of 0.1 /jm 
diameter and 100 /im deep, an aspect ratio of > 1000:1 is 
possible. Applying and addressing electrodes to an array 
of FE nanotubes could generate 3D FRAM structures 
offering high storage density with improved read/write 
characteristics compared to conventional planar stacks. 
On removal of the Si walls, the piezoelectric response 
(expansion/contraction under an applied field) of such 
an array of nanotubes could be utilized for a number 
of MEMs applications. These could include: (1) ink-jet 
printing - delivery of sub-picolitre droplets for lithogra- 
phy free printing of submicron circuits; (2) biomedical 
applications - nanosyringes, inert drug delivery implants; 
(3) micropositioners or movement sensors. 

Almost no theoretical work has been published on the 
physics of ferroelectric nanotubes. Analytical solutions 
for the effects on the dij piezoelectric coefficients of hol- 
low tubes have been given for bo th the case in which po- 
lariza tion P is along the length z ([Ei3eiie^ex_aiigLR ainesb|, 
12003b) and for P radial ijEbenezer and Abrahaml 12002ft . 
they did not however solve the azimuthal case where po- 
larization goes around the tube. It i t this latt e r case 
which been measured as hysteresis bv iLuo et al\ (12003ft 
with a tube lying on a bottom electrode with a semicir- 
cular sputtered top electrode. Important matters such as 
the dependence of Tc upon tube diameter have also not 
been examined. 



VII. CONCLUSIONS 

In this review we have sought to cover the important 
advances in recent years in the physics of thin film ferro- 
electric oxides. At the present point of time ferroelectric 
thin films memory devices have reached a point of ma- 
turity where they are beginning to appear in real com- 
mercial devices. At the same time new directions such as 
the drive to faster, smaller, nanoscale devices and non- 
planar geometries are evolving and new levels of physical 
understanding will be required. Over the next years it is 
expected that first principles computational approaches 
will continue to develop, suggesting a new synergy be- 
tween the computational modeling and experimental re- 
alizations of ferroelectric systems with new and exciting 
properties. 
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